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ABSTRACT 


The results of a program to study the analysis and design of composite 
materials and structures is reported. Emphasis was placed upon three 
major areas: The definition of design criteria for laminates including 

studies of basic failure mechanisms; the definition of unique design concepts 
to enhance the beneficial characteristics of composite materials and to utilize 
them in structures; and the analysis of composite materials property test 
technique s. 


FOREWORD 


This document is the annual report on the program entitled "Study of 
the Relationship of Properties of Composite Materials to Properties of 
Their Constituents". The program was performed for the National Aero- 
nautics and Space Administration under Contract NASw-1377 and was 
monitored by Dr. R. W. Leonard of the NASA Structures Research Divi- 
sion. 
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INTRODUCTION 


Development of methods of analysis of the strength and stiffness of fib- 
rous composites has advanced to the point where it is feasible to establish 
rational initial design procedures for composite structures. These procedures 
are naturally subject to revision as continuing studies enhance the understanding 
of composite failure mechanics. In conjunction with the theoretical develop- 
ments, there is the unsatisfied need for definition of suitable techniques for 
material property measurements. Accordingly, the program described here- 
in treated the above problem areas and the present report presents the results 
of investigations leading to: The definition of design criteria; the enhancement 

of the understanding of failure mechanisms; the definition of unique design con- 
cepts for composites; and the development of improved techniques for com- 
posite property measurement. 

The section on "Design Criteria" describes the failure criteria for lami- 
nates, which have been computerized for structural efficiency analysis. Fur- 
ther studies of the laminae failure mechanisms upon which this is based are also 
described, including some treatment of time dependent behavior. The applica- 
tion of these criteria is treated in the "Design Concepts" section, wherein the 
utilization of the high uniaxial compressive strength of fibrous composites is 
emphasized. Columns, plates and panels of combined composite and metallic 
construction are designed and their potential is assessed. Concepts for im- 
proved material performance including three phase and isotropic three dimen- 
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sional materials are also defined. In the final section, "Studies of Materials 
Properties Tests" the results of the analysis of several current test tech- 
niques are presented along with suggested new methods for uniaxial strength 
measurements. 
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DESIGN CRITERIA 


The parametric evaluations of fibrous composites for aerospace structures 
(e 0 g. Refs. 1 and 2) have indicated the attractive potential of composite 
structures which are configured to achieve high stresses. When efficient 
stiffening arrangements permit the use of the high modulus composite materials 1 
high strength, the resulting structures are shown to be substantially lighter 
than metallic structures. These studies emphasized the need for better 
definition of composite strength. The prediction of laminate strength from a 
knowledge of constituent properties is a complex undertaking and, in general, 
strength estimates most suitable for design are obtained experimentally. 
However, an analytical estimate is required in the assessment of the potential 
of candidate composites which have not yet been brought to the practical fab- 
rication stage. Indeed, the concept of analytically predicting composite pro- 
perties for use in a structural application analysis, is an essential part of the 
search for guidelines for the development of improved composite materials. 

The strengths of fibrous composites, perhaps to an even greater degree 
than their stiffnesses, are complex functions of the anisotropies associated 
with the uni-directional character of filaments. For filaments in one (the 
loaded) direction only, strength intension (Ref. 3), compression (e. g. Ref. 4) 
and shear (e. g. Refs. 5, 6) have been related to the properties of the consti- 
tuents. Further study of these problems is described subsequently. Exten- 
sions to reinforcements in other directions have now been incorporated into 
the computer programs for structural properties of anisotropic composite 
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elements (SPACE). The aim is to provide a strength assessment for prelimi- 
nary design studies. 

Numerous strength theories have been proposed; the validity of which can 
only be justified by experiments which are themselves complex and formidable. 
A comprehensive tabulation of strength theories has been presented in 
Reference 7. Among the theories proposed is the strength theory in References 
8, 9. In this theory, the strength of a fiber- reinforced composite, considered 
as quasi-homogeneous and anisotropic, is governed by a continuous failure 
surface of Hill (Ref. 10). Once the failure surface is determined, the strength 
of the composite body under any type of surface loading can be determined in 
a straight-forward manner. This appears to be a reasonable approach for 
composites with elastic perfectly-plastic fibers and matrix. 

For most composites there is a vast difference between the strength of 
the filaments and the strength of the binder, and failures in the binder may 
be encountered as the maximum stress direction varies from the filament 
direction. In the simplest case, for example, of a unidirectionally reinforced 
composite in tension, if the angle between the tensile load and the reinforce- 
ment direction is increased gradually from zero to ninety degrees, three 
primary failure modes can be expected to be encountered: first, at small 
angles, tensile failure of the filaments; second, at intermediate angles, 
shear failure in the binder; third, as the filaments become oriented mostly 
transverse to the load, tensile failure in the binder. These failure modes are 
essentially independent of one another. Particularly for the change from 


4 


tensile failure in the filaments to shear in the binder, there appears to be 
little reason to expect a gradual transition of the type that leads to a smooth, 
"yield surface" for homogeneous materials. 

When the various failure modes are independent of each other for a com- 
posite having an oriented structure, the applicability of a continuous function, 
like Hill’s anisotropic yield condition for a homogeneous material (Ref. 10), 
to represent a yield or strength criterion appears open to question. Accordingly, 
the approach utilized herein to determine strength criteria for composite 
laminates has been to determine separately the strengths for all possible 
failure modes. Thus, to a degree, a family of failure surfaces representative 
of the material will be generated, and the lowest of them for any loading 
condition will be the governing one. 

This approach is described in the following sub- section. This is followed 
by a presentation of the recent studies of the strength properties of a uniaxial 
laminate. Note that these latter quantities are required both for the dis- 
continuous and the continuous failure surface models. The analysis is developed 
first for composites where both fiber and matrix are elastic-plastic and a 
continuous failure surface is defined. These results are then specialized for 
rigid brittle fibers to generate strength values which can be used in a maximum 
stress failure theory. 

Laminate Strength 

The strength analysis of a laminate of layers of uniaxial fibrous compo- 
sites utilizes the elastic analysis, under given surface loadings, of the state 
of stress in each laminate layer considered as quasi-homogeneous, i. e. 
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locally heterogeneous, but grossly homogeneous. If the surface loadings 
increase monotonically and proportionally, there will be a stage at which 
the stress in one (or more) layers of the laminate is at a failure point and 
the layer, being assumed to fail, is replaced by a new degraded layer having 
an assumed mode of degradation. As successive failure of constituting 
layers proceeds, a redistribution of stress among the laminae occurs and 
the slope of the load-deflection curve is discontinuous. The ultimate strength 
of the laminate is reached when all the constituting layers have failed. 

In each layer, the stresses referred to the principal axes of anistropy 
are computed. If the shear or transverse stress is equal to its corresp ending 
yield stress, the lamina is considered to hold that stress level for those com- 
ponents and to have additional stiffness only in the fiber direction. As 
successive failure ot constituting layers occurs, the entire load-deflection 
history can be traced until a failure in the fiber direction occurs or until all 
layers yield, at which point the associated applied load is defined as the failure 
load. An illustration of the application of this approach is presented in Fig. 

1 where calculated stre ss- strain curves for two simple laminates of E-glass 
in epoxy are shown. For more general laminate configurations there will be a 
greater number of straight line sections in the stre ss- strain curve. 

For those applications where only the limiting stress levels - and not the 
entire stress-strain curve - are required, a simpler approach to the defini- 
tion of failure criteria appears reasonable. This approach is based upon the 
concept that the first departure from elastic behavior is a most significant 


6 



point on a composite material stress- strain curve. It is desired to keep 
actual stress levels below this point in a fashion analagous to the use of the 
yield stress for metal construction. Similarly, at the ultimate stress level, 
the transverse properties of the individual layers have generally deteriorated. 
Thus this level can be approximated by using a "netting" analysis with the 
uniaxial strength properties of the individual layers. The application of these 
principles in the definition of failure criteria is described below. These 
methods have been incorporated as a subroutine of the SPACE computer 
program. 

The basic stress strain relations for the laminate are given by: 



Where N. and M. are the laminate stress and moment resultants, e are 
11 oi 

the middle surface strain components, R. are the curvatures, and a , b 

1 ij ij 

and d are the laminate elastic constants obtained by suitable integration of 
the laminae elastic constants. For definition of laminate strength we desire 
to treat the stresses arising from a set of applied stress resultants when the 
curvatures are prevented. Thus we can consider above equation rewritten as: 

N 


/ _ V 

r -1 ' -1, 

( 6 ) 

a j -a b 

} ° 

i 

i 

} 

J [ = 

\ 

i 

1 \ 

-i ; -i 

V m ; 

_ ba d-ba 
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is the inverse of the matrix a. 


where a 


-1 


For zero curvatures we have 


{e} = La ] _1 {n} 


and 

- ( k \ - ( k ) - 

or. = C. e 
1 ij oj 

where O' are the stress components in the k layer referred to the 
laminate axes, 

(k) 

C ^ are the elastic moduli of the same layer referred to the same axes. 

From these equations we may find the stresses, crT\ referred to the lamina 

1 

principal axes from: 


a' k) = T..O k )i (k > 

1 IJ J 


no sum on k 


where T is the transformation matrix for rotation of coordinate axes. 

ftk . 

P is the angle from laminate principal axes to lamina principal axes. 


_ i k 

To evaluate laminate strength, the stress components CJ are evaluated for a 

laminate having its laminae thicknesses and the total load normalized. Thus, 

the total thickness is unity and the load vector for axial load, for example, is 

(1, 0, 0). With the stress components known, and the maximum strengths, 

(k) 

a iu defined for each of the stress components, (as in Refs. 3, 4, and 5) 

the ratios (k) 

0 ' ' 
iu 


(k) 

a 

iu 
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are formed. The lowest ratio of this set of 3n quantities (for an n layer 
laminate) is the material yield stress, CT . 

y 

The same analysis is now repeated for the case where the moduli E^ , 

^ 12 7 an< ^ ^21 are se * ec l ua ^ *° zero for all layers. Thus only E^, the modulus 

in the fiber direction is non-zero and we have (for E^ = v^E^) a ''netting 1 ' 

analysis. Here there are only n stress ratios and the lowest one is taken as 

the material ultimate stress, cr . 

u 

As an application, the yield stress of a symmetric biaxial laminate sub- 
jected to an axial tensile load was treated. The results for these laminates 
are presented in Fig. 2. Each curve represents the results generated for one 
of the three stress components and the lower envelope curve is the design 
yield stress curve. 

Lamina Strength 

It has been postulated (Refs. 8 and 9) that the strength of a unidirectional 
fibrous composite, considered as quasi-homogeneous and anisotropic, can be 
represented by a surface having the form of Hill's Generalized von Mises' 

Yield Condition (Ref. 4), namely: 

2£ ( T iP ■ F ( T Z 2 - T 3 3 ) 2+ G ( T 33- T ll) 2 + H ( T 11- T 22 ) 2 + 2LT 2 3 <» 

+ 2M T 2 +2NT 2 =1 

31 12 

where the coefficients F , G, H, L, M, N are parameters characterizing the 
state of anisotropy and T are components of the stress tensor referred to the 
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principal axes of anisotropy x 1 x 2 x 3 where the Xj axis is parallel to the 
fibers.* In general, the yield condition ( 1) can be represented by a surface in 
a stress space. A stress point within the yield surface represents a "safe" 
stress state. Yielding can occur if the stress point is on the yield surface. 
Since unidirectional fibrous composite layers in laminates are thin compared 
with their lateral dimensions, only T n> T ^ and T are considered non-zero 
if the x^ axis is along the thickness of the layer. Furthermore, since uni- 
directional fibers are randomly located in a composite layer, it can be assumed 
to be transversely isotropic. Then the yield condition (1) reduces to 



where and X^ 2 are the normal yield stress in the direction of the 

fibers, the normal yield stress in the direction transverse to the fibers and 
the yield stress in axial shear of the composite, respectively. These are the 
three basic strength characteristics of the unidirectional fibrous composites. 
Once these are known, the yield condition (2) can be employed to determine 
whether a combined state of (plane) stress can cause failure of the composite. 

In what follows, effort is made to evaluate the quantities, Xj , X 2 and 
^12 analytically in terms of the strength and geometry of the constituents. 


* Henceforth, unless otherwise specified, i, j = 1, 2, 3; Summation on 
repeated indices is implied. 
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Limit Analysis of Unidirectional Fibrous Composites 

The composite material under consideration consists of a 
relatively soft matrix material in which stiffer fibers are embedded. Initially 
both materials are assumed to be elastic-perfectly plastic and satisfy the von 
Mises 1 yield criterion. As shown in Figure 1, referred to an orthogonal 
Cartesian co-ordinate system, a typical unidirectional fibrous composite is 
taken to be a cylinder with rectangular cross-section. Circular fibers run- 
ning from base to base of the specimen are in x -dire ction. Limit analyses 
of such a specimen with various arrangements of both elastic brittle and elastic 
perfectly plastic fibers in an elastic-plastic matrix under various types of 
surface loading will be described. 

Elastic-plastic Constituents 

In this study only the "random array" geometry is considered: circular 

fibers of various diameters are randomly located in the specimen. Each of 

them can be surrounded entirely by a concentric cylinder of matrix material. 

A cylinder consisting of a fiber of radius r^ and the outer matrix- shell of 

radius r^ is called a composite cylinder. It is assumed that a constant 
r 

g _ can be chosen so that the composite cylinders are non-overlapping. 

r b 

The entire specimen can then be considered as an assemblage of composite 
cylinders and the ramining matrix volume. The lateral boundary of the spe- 
cimen may touch or cut through some fibers. In both cases the associated 
composite cylinders are "incomplete 1 ' . Since in practice fiber diameters are 
very small compared with the transverse dimensions of the specimen, the 
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total number of fibers in the interior of the specimen is much larger than the 
total number of those possibly on the lateral boundary. Hence, the total 
volume of "incomplete" composite cylinders is much smaller than that of the 
"complete" cylinders. 

If Vj and V 2 denote, respectively, the total volumes of the specimen, 
the composite cylinders and the remaining matrix in the specimen, the 
following obvious relation holds: 


V = V 


1 


+ V. 


(3) 


In the case where the entire specimen is occupied by composite cylinders, 
V 2 = 0# following Hashin and Rosen (Ref. 11) this distribution is called the 
"random array". 

The von Mises 1 yield criterion which the fiber and matrix materials are 
assumed to obey has the following form (Ref. 12): 


S.. S.. 
2 


(4) 


where S are components of the stress deviator and k is the yield stress in 

simple shear for the fiber material (denoted by k^.) or for the matrix material 

(denoted by k. ). 

b 

Under the conditions of plane strain perpendicular to the x^-axis, von 
Mises 1 yield criterion (4) reduces to 


r. _ 

22 


33 


4T. 2 

Li 


4 k 


(5) 
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where and are components of the stress tensor in the transverse 

plane. 

In order to evaluate the three basic strength characteristics (Refs. 8, 9) 
of the composite specimen shown in Figure 3, the specimen is subjected to 
axial shear stress transverse tensile stress and longitudinal tensile 

stress T| j , respectively. The upper and lower bound theorems of limit 
analysis of plasticity (Refs. 13, 14) will be used to obtain bounds for the limit 
loads T 12 , T ^ and T ^ which represents the lamina strengths. 


Case 1 . Axial shear stresses T applied on the boundary of the composite 

specimen. 

According to the lower bound theorem, a uniform shear stress field T 
can be chosen as the statically admissible stress field. Since can nowhere 

violate the yield condition (4) for both matrix and fibers, it follows that a 
lower bound for the limit load T ^ is 

T 12 )l = k b 

For upper bound construction, a kinematically admissible velocity field 
is chosen as follows: 

(a) In (and thus also on the boundary of the composite cylinders) and 
on the entire lateral boundary of the specimen, 
u. 




b x i 


(6) 


u 3 ) ( 0 ) 

where is a positive number. 
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(b) In any composite cylinder, referred to a local cylindrical 


polar coordinate system, * 


u. 


- Y r cos 9 


u 


u c. . 
\ ^ / 


Yj x^ cos 9 

- Y. x, sin 6 
i i 


for 0 < r < r 


f 


V (1 V 

? \ 1 - — 7- r cos o 

1 -r r 2 


Yj x cos 8 


for r < r < r 
f b 


- Y jXj sin 


* The velocity field in any "incomplete" composite cylinder is defined by 
solving similar elastic displacement boundary value problems for the 
"incomplete" composite cylinder. However, since the volume of the "in- 
complete" composite cylinder is small, the difference between their actual 
contribution to the dissipation function and that obtained by treating all 
cylinders as "complete" is negligible. This approximation is implied in 
the subsequent analysis wherever a similar situation arises. 
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The velocity field (7) in a composite cylinder is the elastic displacement 
solution to the boundary value problem with the boundary condition (6) pre- 
scribed. The problem is the same as that which was formulated by Hashin 
and Rosen (Ref. 11) with the modification that for 0 < r < r , the velocity field 
is associated with rigid body motion. 

With this velocity field constructed for the entire specimen, the dissipa- 
tion density function and the rate of external work done can be obtained to 
yield an upper bound for T ^ (Ref. 5). 


12 / U 


1 2 TT 


f | 


p 4 \ e 2 

R / I 1+ — — J+ 2 — — cos 0 d 6 dR (8) 


(. „2\ 

V U-P ) J J \j \ R" / R 
P o 


The above expression is for "random array" in which = 0 and the fiber 
„ 2 

volume fraction v^ = |3 . 

The integral in (7) is integrated numerically for different fiber-volume 

(t L \ 

fractions. The result is shown in Figure 4 where \ 12 ) v in (8) is plotted 

2 \ 

as a function of v f (0 < v f < 1, v = £ ). Note that particularly in Figure 4, 


lim 
v^ - * 0 


lim 
V f ' 


( T 12 h 


1. 


4_ 

IT 


F rom the above result, it is concluded that is at most about 27% above 

and at least the same as the yield stress in shear for the matrix. 
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Case 2 . Transverse tensile str 


esses T 22 applied on the boundary of the 


composite specimen. 

For lower bound construction, a uniform tensile stress T ^ throughout the 
specimen is chosen as a statically admissible stress field. Since the von 
Mises' yield condition (4) can nowhere be violated in the specimen, the lower 
bound associated with the constructed statically admissible stress field j 8 

For upper bound construction, a kinematically admissible velocity field is 
chosen as follows: 

(a) In V 2 (and thus also on the boundary of the composite cylinders) and 
on the boundary of the specimen. 



where Y 2 is any positive number. 
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(b) In any composite cylinder, referred to the local coordinate system, 



2 x. 


Y, 

- 2 x 


2 2 
A l +A 2 * 2 + 3 X 3 


2r 


f 


+ A, r f ( - x 2 + 3x 3 > + A r f (x 2 ' x 3 )_ 

4 — . 

2 2 


o 2 2 4 2 2 2 _ „ 
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The velocity field (11) in a composite cylinder is the incompressible 
elastic displacement solution to the boundary value problem with the boundary 
condition (10) prescribed. The problem was formulated by Hashin and Rosen 
(Ref. 11) in evaluating the effective plane strain shear modulus of elastic 
fiber-reinforced composites. The constant p in (11) considered as the displace- 
ment solution to the elastic problem is the ratio of the elastic shear moduli 
of the fiber and matrix materials. However, it can be considered as merely 
a parameter when (11) is used as a kinematically admissible velocity. 

For the case of the "random array", an application of the upper theorem 
gives an upper bound 



( 12 ) 


where 
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The right hand side of expression (13) was obtained in Ref. 5 for 

an upper bound of the limit load (normalized with respect to k^) for transverse 

shear stresses applied on composites reinforced with rigid fibers. In Figure 

(r L \ 

5, \ 22 / v is plotted as a function of v according to (13). It is a monotonically 

2k b 


Jd l l 

increasing function of v. with lim [T__ ) TT . . lim \ 

f v _ o V 22 L U = 1 and [22 [ U = » . 

f 2 h, f 2 k b 

On the other hand, for any finite 7], (t^ ^ ^ ^ can be plotted as a function of 


:1 % 


(r,h, 


by evaluating numerically 1^ and 1^ in (12). For example, for r\= 100, x 22 7 U 

2k b 
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is also a monotonically increasing function of v 


lim 


L\ 

T ”) 


22 ) U - 


( T ^ L )u \ 

£ 2k b ' k b ) ' 


lim 
v 


2k b 


which is shown in Figure 3, with f = 10 . 

k~ 


and 
It is 


observed that for fiber volume fractions smaller than about 75%, the upper 
bound for the limit load T ^ is higher that for T)- « but for higher fiber volume 
fraction, the opposite is true. Therefore, for any fiber volume fraction between 
0 and 1, the lower value of the bounds obtained from r\= 100 and T]-** will give a 
better estimate of and it is shown in Figure 5 in solid line as a function of 


It is interesting to note that for r\ = 1, (12) reduces to the following 
simple form 



(14) 


which is commonly known as the "rule of mixtures". From Figure 3, it is 

seen that the straight line represented by (14) will be higher than the chosen 

curve for ( 22 )u for 0 <v -=j. Therefore, it is concluded that the "rule of 

2k b 

mixtures to be used in this case would overestimate the composite strength 

(r L ) 

In Figure 5 , \2Z / L obtained from (9) is also shown. The difference 
2k b 

between the upper and lower bounds for v -0 is due to the fact that in the 
upper bound construction, a plane strain velocity field is used as a kinemati- 
cally admissible velocity field. 

C ase Longitudinal tensile stresses applied on the boundary of the 

composite specimen. 
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The study of longitudinal strength of unidirectional fibrous composites is 
extensive. Various models and failure mechanisms have been proposed in 
the literature. Here, the yield strength in axial tension is obtained by the 
construction of very simple velocity and stress fields, for composites having 
elastic-plastic fibers of uniform strength. 

For upper bound construction, in the entire region of the composite spe- 
cimen, a kinematically admissible velocity field is chosen as 



An application of the upper bound theorem gives 


( T n) = T f 


v, + T v 
f b v 


(15) 


or 


tli)= I +v/— - l\ 


where T and T are the yield stresses in tension for the fiber and matrix, 
t b 

respectively and v^ = 1 - v^. 

Relation (15) is known as the "rule of mixtures" for the predication of the 

(r L ) 

strength of a composite. If 7^ and are given, \ 1 1 / U is a linear function of 

T b 

fiber-volume fraction v^ . 
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F or lower bound construction, if the applied stress r is assumed to be 

11 

uniformly distributed on the boundary of the composite specimen, then a 
lower bound ^ ^ L is equal to T b . However, if we assume that the tensile 

stress T u is not uniformly distributed on the boundary surface and we are 
only interested in the average stress intensity on the boundary surface that 
causes failure, then a higher lower bound can be obtained. In fact, a statically 
admissible stress field can be chosen as follows: 

In the region occupied by fibers 


T. 

ij 


T 0 0 

0 0 0 


0 0 0 


(16) 


In the region occupied by the matrix. 


T .. 
ij 


T b 0 0 

0 0 0 

_0 0 0 


(17) 


It is obvious that the stress field expressed in (16) and (17) satisfy the re 
quirements to be statically admissible. 

The average traction corresponding to this stress field is therefore 


Vf + r b v b 


(18) 


which can be taken as • 


According to (15) and (18), the upper and lower bounds for 7 L 


11 


coinciue. 


Therefore, this strength can be determined from the "rule of mixtures", for 
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this type of composite, (i. e. elastic-plastic fibers and matrix of uniform 
yield strength),. 

From the above results, it is observed that obunds for are far apart, 
especially for high volume fractions. Further effort should be made to im- 
prove the bounds in order to have a better estimate of the limit load. For lower 
bound construction, uniform stress distribution used as statically admissible 
stress fields can only give lower bounds which are the corresponding matrix 
yield stresses. In order to obtain higher lower bounds, one has to assume 
applied tractions to be non- uniformly distributed on the boundary surface 
according to the properties of the fibers and matrix,, Equilibrium stress fields 
can then be constructed in equilibrium with the applied tractions. Then the 
lower bound theorem can be applied to obtain higher lower bounds. To decide 
the distribution of the applied tractions on the boundary and to construct an 
equilibrium stress field in the body is not at all easy in general. The success 
in the construction of such a statically admissible stress field for Case 3 is 
due to the simplicity of geometry and loading conditions. 

Brittle Fibers 

In the previous section, the strength of unidirectional fibrous composites 
has been evaluated by obtaining bounds for the basic strength characteristics. 
The theory assumes the existence of a continuous failure surface. This hypo- 
thesis appears reasonable for composites with elastic-perfectly plastic 
fibers and matrix such as metal fibers and matrix. However, for contemporary 
high strength fibers embedded in epoxy resin matrix, this n elastic-perf ectly 
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plastic model" may seem inadequate. Instead, it seems more suitable to 
assume elastic-brittle fibers and an elastic-perfectly plastic matrix which 
obeys the von Mises' yield criterion. Upper and lower bounds have been ob- 
tained for the following types of surface loadings: 

Tractions equivalent to a uniform shear stress T applied on the en- 

1 W 

tire boundary surface of the specimen. 

2. Tractions equivalent to a uniform shear stress T applied on the 

<-i J 

entire boundary surface under the conditions of plane strain. 

3. Tractions equivalent to uniform uniaxial tension T applied on the 
entire boundary surface under the conditions of plane strain. 

4* Tractions equivalent to biaxial uniform tension T and T 

22 3 3 

( T 22 ^ T 33^ a PPli e d on the entire boundary surface under the conditions 
of plane strain. 

5. Tractions equivalent to combined in-plane shear T and transverse 
tension {T = OL >0 ) applied on the entire boundary surface. 

Since the method of analysis is similar to what has been presented in the 
preceeding pages in obtaining bounds for the basic strength characteristics, the 
details of analysis will be presented in Appendix A. However, results are 
summarized as follows: 
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(r L ) 

for the "random array" geometry. The above expression for \ 12 /U is the 

same as ( 8 ) which is shown in Figure 4 as a function of (0 < ^ <1, = £ ) 

since the same velocity is used here. 

For arbitrary geometry of arrangement of fibers, (^12 ) U _ ] + (~ \\ 

L x - ' v ' 

(t \ 

which is higher than \ 12 / U given by ( 8 ). 
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for the "random array" geometry. The above expression for ^7^ ^ ^ ^ is the 

k b 

same as (13) which is plotted in Figure 5 as a function of with 


‘r-o 

f K 


and 


lim 

vF 


i ( Tz3 ) u 


since the velocity field chosen here can be obtained from ( 10 ) and ( 11 ) through 
an orthogonal transformation together with the limiting process of r\ co 
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for a >0, 


For the case of the "random array” geometry, 




\ + ““ 


and 





(19) 


( 20 ) 


26 


whe r e 


1 27T 
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3 p R 


2 2 (4P 4 + ^+i) 


Numerical calculation is performed to obtain and (r from 

^ 1 22 /U V 12 ;u 


(20) for different values of P and & . The results are summarized in Figure 6 

2 

in which f3 = 0. 8 is the highest fiber volume fraction shown. The dotted line 

represents (19) which gives the lower bound for (tv ^and (j ^\for any p . 

\ 22 / \ 12 / 

It is worth mentioning that for this model of elastic -brittle fibers embedded 
in elastic -perfectly plastic matrix, in the upper bound expressions, only the 
matrix strength and fiber -volume fraction appear - the brittle strength of the 
fibers is not involved. This is due to the fact that rigid body motion of fibers is 
always incorporated into the kinematically admissible fields. Hence, regions 
occupied by fibers in the composite body contribute nothing to the dissipation 
function. The result is reasonable since in reality, elastic moduli and brittle 
strength of fibers are much higher than the corresponding elastic moduli and 
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strength of the matrix material. On the other hand, the longitudinal strength 
along fiber-direction of fibrous composites under tension or compression cannot 
be obtained by limit analysis of plasticity because of the presence of the high 
modulus, high strength elastic -brittle fibers. Instead, a statistical failure 
theory was established for tensile strength [3] and a fiber -buckling model was 
postulated for compressive strength of fibrous composites [4] and [l 5 ]. These failure 
mechanisms are based on experimental observation. Furthermore, as men- 
tioned previously, for this mathematical model of composite material, the 
existence of a smooth yield surface that governs the strength of the material is 
still an open question. Based on experimental investigation, the failure mech- 
anisms for individual cases seem to be different and independent of one another. 
Accordingly, the approach to the determination of the strength criterion for uni- 
directional fiber-reinforced composites as well as fibrous composite laminates 
is as described earlier. Instead of using a continuous yield surface, all possible 
failure modes are considered. To determine which one is dominant for a parti- 
cular type of loading condition can sometimes become a tedious task especially 
for fibrous laminates under various types of design loading. However, this 
approach has the advantage that the entire load -deflection history can be traced 
until complete failure occurs. With modern high speed digital computers avail- 
able, a systematic strength analysis can be programmed to obtain accurate 
results . 

For certain composites, it is possible to utilize, as an alternate approach, 
the application of limit analysis methods to the laminate as a whole. This is 
treated in the following section. 
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Application to Laminates 

It is interesting to note that the kinematically admissible velocity fields 


constructed to evaluate the upper bound for the three basic strength characteris- 
tics can be used to obtain upper bounds of limit loads for in -plane loading applied 
on laminates . 

To demonstrate the method, consider a laminate composed of uni- 
directional fibrous laminae subjected to in-plane shear stress j . The 

1 Ca 

L 

problem is to find an upper bound for the limit load T . A typical constituting 

12 

layer, the k^ k layer, in which fibers are all running in x^ direction, is shown 

in Figure 7. Depending on the orientation of fibers in the layer, the principal 

axes x x x can be defined by an angle 0^ k ^ measured from the lami- 

12 3 

nate axes x x x . 

1 2 3 

For upper bound construction, a kinematically admissible velocity field 

is constructed in the following manner: 

The velocity field (p , p , p ) = ( 0, V x 0) is assigned in V 

12 3 11 4 

of each layer , referred to x ^ x^ x^ axes . In the k^ k layer ref erred to its prin- 
cipal axes x x x„ , the above velocity field is transformed into the 
* 12 3 ’ 


following form: 

- (k) 

Hi 

P ,k) 


v .• fi< k > n< k > < k > _ v .2 a< k > 00 

Y sin 0 cos □ x - Y sin d x 

X 11 u 

v 2 (k) (k) . fi (k) ,(k) (k) 

Y cos x - Y sin b cos b x_ 

111 L* 


( 21 ) 


- (k) 

P. 


The right hand side of equation (21) can be decomposed into four parts 


so that 
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( 22 ) 



where 





and 



i; (k) + M (k| +-n (k| +t; < k) 

->1 -»2 _>3 A 


0 
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sin 6 x 
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0 


Except for multiplicate factors and reference coordinate system, 
- (k) - (k) - (k) 

u j » u o > and u are the same as those velocity fields constructed 

-» -» 

in V 2 f ° r the three different cases of loading in the proceding section. 
Furthermore, 
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-4 W - 


v • 2 fl( k > W 

-Yj sin b 
0 
0 


0 

• 2 0 (k) (k) 

-Y sin 0 ' x 


V • 2 fl( k > < k > 

-Y j sin b x^ 

.2 (k) (k) 

+ \ Yj sm b X j 


0 


(23) 


— (k) 

where the first part of the decomposition has the same form as u while the 

-► * 

second part represents a rigid body motion of V as a whole. 

The velocity field chosen in is also composed of four parts: The first 

— (k) 

part is due to u , which is the same as (7) with the modification that Y 

2 (k) — (k) 

cos 9 . The second part is due to u^ which is the same as (11) with the 


(k) - (k) 

modification that Y is now replaced by -Y sin G v 1 cos 9 V . The third part 

c* 1 


i(k) 


is due to which is exactly the same as u in V « The fourth part is due 

^3 3 2 

— (k) 

to u^ which can be obtained in the same way as in Case 1 in the preceding sec- 
tion. However, a moment's reflection reveals that the strain rate associated 

— 2 „ (k) 2 (k) 

with u 9 if Y 1 cos 9 V 'is replaced by -Y, sin 9 1 since the additional rigid 

body motion contributes nothing to the value of the strain rate. Now, after 

the kinematically admissible velocity field is constructed as described above, 

the associated strain rate and dissipation density function can be calculated 

without difficulty. Moreover, the rate of external work done can also be cal- 

L 

culated so that an upper bound for T ^ can be obtained. 

As an example, consider a laminate composed of two unidirectional 

(1) ~ (2) _ tt 

fibrous composite layers of equal thickness with 9 =0 and U - . 

In this simple geometry of lamination, the kinematically admissible 
velocity field constructed above reduces to a very simple form. 

For layer 1, the velocity field (21) reduces to 
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For layer 2, the velocity field (21) reduces to 

(25) 

Then, following the principles described above, it is easy to obtain the velocity 
fields and the associated strain rate fields for both layers. After some mani- 
pulation, it turns out that 



(T 


L 


12 


1 

TT(1-P 2 ) 



R 



3 

4 



cos 6 d0 dR 


for the "random array" geometry. 

The above equation is exactly (8), which means that the upper bound 
T 1 2 U 1S n0t hl § her for this l yp e of cross-ply laminates than for the uni- 
directional fibrous composite. 

Similarly, it is easy to show that the three types of velocity fields con- 
structed in the preceding section can be used to construct upper bounds for 
limit loads for laminates subjected to any in-plane stresses. 
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Tensile Strength 

The high-strength, high modulus fibers which are of interest for 
use in composite materials are generally brittle, having tensile strengths 
that must be characterized statistically. Any theory for the tensile strength 
of composites containing such fibers must take into account the dispersion in 

their failure stress levels in order to have any relevance. 

( 16 ) 

Parratt noted the dispersion in fiber strength and suggested that 
failure of a fibrous composite subjected to tensile load occurs when the 
fibers have broken up into lengths so short that any increase in applied load 
cannot be transmitted to the fibers because the limit of interface or matrix 
shear has been reached. 

A theory has been presented (3) for the failure stress of composites 
containing continuous, uniaxially- oriented, brittle fibers in a ductile matrix 
loaded parallel to the fiber direction. This theory predicts that due to the 
distribution of flaws or imperfections in the fibers there occurs a series of 
randomly-distributed fiber fractures as the applied stress level is increased. 

It is argued that a portion of the broken fiber in the vicinity of the 
fracture is ineffective in resisting the applied load. Assuming that the stress 
in a broken fiber is uniformly distributed among the unbroken fibers in the 
cross-section and that this overstress acts over a length equal to the "in- 
effective" length, it is predicted that failure occurs when a weakened cross- 
section cannot sustain an increase in load. In effect, the theory pre- 
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diets that the composite has a strength equal to that of a bundle of fibers 
whose length is the "ineff e ctive 1 ' length. 


The present study considers fibers having a statistical strength dis- 
tribution resulting in fractures at various stress levels as the applied load 
increases. It is assumed that the stress in the broken fiber is distributed 
unevenly to me omei lioci 5 in me cross-section which has a length equal to 
the ineffective length. As a first approximation the effect of this overstress 
is presumed to affect only those fibers adjacent to a break. Failure is as- 
sumed to occur due to an increasing probability of fracture in the fibers ad- 
jacent to a prior break. 


Description of the Model 

The model consists of a two-dimensional composite of length L 
consisting of a ductile matrix in which are imbedded N continuous brittle 
fibers whose orientation is parallel to the applied tensile load. The fibers 
are considered to be composed of layers of length 5 . The total number of 
layers being M = L/6, (See Figure 8. ) 

The quantity 6 represents some length over which the stress is per- 
turbed in the area of a fracture. It is variously referred to as the ineffec- 
tive length or twice the transfer length and several formulae have been proposed 
for its evaluation. Two of the definitions ( 31, (17)are based on an elastic shear- 
lag type of analysis. The axial stress in a broken fiber is found to be 


34 



( 1 ) 



where 


2 


T 






E^ = Young* s modulus of fiber 

= Shear modulus of binder (matrix) 


r 

f 

v 

f 

X 


Fiber radius 

Volume fraction of fibers 
Distrance from end of broken fiber 


a 

0 

fiber end. 


Extensional stress in the fiber at a large distance from the 


It should be noted that in Reference 3 the factor of 2 in the expres- 
sion for T] was incorrectly omitted. The first author defines 6 as the value 

of x for which the stress in the fiber has reached 90% of a , the stress at a 

0 

long distance from the fiber break. On the other hand, Friedman defines the 
ineffective length by means of an approximate step -function stress distri- 
bution which has the same average stress as the distribution of Equation 1. 
Furthermore, this author includes the ineffective portion on both sides of 
the break whereas the first definition uses only one side. 

The expressions for the two ineffective lengths discussed above 
are, respectively 
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f 



1-v 


1 / 2 ' 


1/2 


1/2 
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F 


d 



1-v 


1/2-1 


1/2 


1/2 


( 2 ) 


( 3 ) 


L f J 

If the shear stress between the matrix and fiber is assumed to be 
constant, as in the case of plastic flow or frictional shear stress due to inter- 
face failure, the fiber tensile stress becomes constant at a finite distance 
from the fiber end. This distance is called the transfer length l by some 
authors and is given by the expression 
CT r 

■L = -°J- 

t 2 r 


where T is the (constant) plastic shear stress of the matrix, or, if there is 
interface failure, the frictional force between the matrix and fiber. 

The ineffective length can also be determined experimentally by a 
photoelastic examination of the stress patterns in the vicinity of a fiber break. 
In this investigation the experimentally determined ineffective length is used 
where available. 

Statistical Analysis of the Model 

It is assumed that the strength of the population of fiber elements of 
length 6 can be characterized by a cumulative distribution function F(a). That 
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is, the probability that an arbitrary element has a failure stress level less 
than, or equal to a is F(a ). 

The expected number of elements in the composite that will fail under 
fiber stress a is 


E = MNF(ct) . (4) 

This expression includes the possibility of further fractures of adjacent 
fibers. If the composite does not fail, the stress in the broken fibers is dis- 
tributed to the other fibers in the cross-section in a complex manner. As a 

/ no \ 

first approximation, Hedgepethr used a shear lag analysis to determine the 
average stresses in fibers adjacent to an arbitrary number of broken fibers. 

He considers an infinite two-dimensional array of fibers subjected to 
tensile load parallel to the fiber direction which is uniform at a great dis- 
tance from the fracture area. The ratio of stress in the two fibers adjacent 
to a run of r broken fibers to the uniform applied stress at infinity is, for a 
static stress distribution 


= 4- 6- 8* * (2r+2) 

r 3* 5° 7* (2r+l ) 


(5) 


Hedgepeth calls K a stre s s - concentration factor, but in this paper it will be 
referred to as an over st res s -factor so that it will not be confused with stress 
concentration factors found by an "exact" analysis of the stress distribution. 

For the case where r fibers break simultaneously the author demon- 
strated that the ratio of the maximum dynamic stress to the static stress in 
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the fibers adjacent to the break increases from 1.15 for r = 1 to a limit of 
1. 27. The values for r = 2 and r = 3 are 1. 19 and 1. 20 respectively. 

In the present paper it is assumed that the overstress in two fibers ad- 
jacent to the broken ones exists over the entire ineffective length. Therefore, 
the probability that an element adjacent to r broken elements will fail is, 

approximately, the probability that its strength lies between a and K a . 

r 

This probability is equal to 

F(K r a) - F(a) (6) 

This approximation is justified, as will be shown later, by the fact that com- 
posite failure occurs for small values of F(a) and therefore the probability of 

having adjacent fractures because both fibers have strength less than C, which 

2 

probability is proportional to F(o) , and the probability of interaction of frac- 
ture groups is relatively small. 

Given that a single element is broken the probability that one of the two 
adjacent fibers will break is 

2 

p 2/1 - 2 [f(K 1 ct ) - F(ct )] - 2 [f(K^ - F(ct )J . (7) 

The probability that both adjacent fibers will break simultaneously is 

P 3/l = [ F(K l o) - F <°>] 2 (8) 

It should be noted that each of these expressions does not exclude the possi- 
bility of further fractures. 
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It is now assumed that only the two fibers immediately adjacent to a break 
are subjected to an overstress and that all of the remaining fibers in the cross- 
section have a stress level equal to the average stress CJ . If one of the fibers 
adjacent to a single fracture breaks the fibers adjacent to the two broken 
fibers are subjected to a stress level K^CJ. One of these overstressed fibers 
was previously exposed to a stress level K^CJ while the other saw only the 
average stress a . The probability that one of the two fibers will break is 


[f(K 2 <7) - F(K o)j + [F(K tr) - F(CT)] 

-2 [f(K a) - F(K d)] [f(K a) “ F < a )] 


The probability that both fibers will break simultaneously is 

P 4/ 2 = [ F ( K 2 CT ) - F < K 1 CT )] [ F < K 2 a) - F < a )] (10) 

If both fibers adjacent to an initial fracture break there will be three 
broken fibers in a row and therefore the two fibers adjacent to this group, 
which were previously at a stress level a, will be subjected to a stress K^g ' 
Again it is possible for one or two of these fibers to break, and so on. 

By now the process and complexity of the problem should be evident so 
that the expressions for further fracture probabilities are presented without 
discussion. Note that that there are two paths by which a state of three 
broken fibers can be reached from a single break; A) by the simultaneous 
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breaking of both fibers adjacent to the initial break, B) by the successive 
breaking of two fibers. 

The expression p.^ represents the probability of having i fibers 
broken given that j are already broken. The letter z represents the particular 

path if there is more than one. The letters A and B refer to the paths des- 
cribed above. 


P 4/3A = 2 [ F < K 3 a ) " F(ct)] - 2 [f(K 3 CT) - F(<j)] 2 

P 5/3A = [ F (K 3 a) - F(a)] 2 (lla-d) 

P 4/3B = [ F(K 3 a) " F(K 2 Q) ] + [ F < K 3 a) ' F (^)] ~ 

2 [^F(K 3 ct) - F (K^a)J |^F(K 3 ct) - F(a)] 

P 5/3B = [ F(K 3 CT) - F < K 2 a >] [ F ( K 3 °) - F (a)] 


Probability of Cumulative Fractures 
Each of the E fracture sites distributed throughout the composite is a 
nucleus for further fiber breaks because of the overstress in adjacent fibers. 
The probability that an element will break followed by the fracture of at least 
one adjacent element is 


F(g) (p. 


+ p ) . 

H 3/r 


( 12 ) 


Therefore, for the composite as a whole the expected number of groups of 
two or more broken fibers is 
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(13) 


M(N-l)p. 


2 " 2 


and the associated probability of having at least one such group is 


p, = Hl-ip ) M(N - 2 » 
2 2 P 2 


(14) 


where the factors of 1/2 are introduced to account for the independence of or- 
der of two fractures. 

The probability that an element will fracture followed by the breaking of 
at least two other fibers in a row is equal to the probability of at least two 
fractures less the probability that two will break without further fractures 
occurring. This probability is 


P, = P, “ F ( a ) P 


2/1 


( 1 - 


P 3/2 _P 4/2 > 


(15) 


The probability of having at least one group with three or more frac- 
tures is 


i n 1 . M( N - 2 ) 

P 3 = 1 - (1 ‘T P 3 


( 16 ) 


and the expected number of such events is 


e 3 = — P 3 M(N-2) 


(17) 


The analogous expressions for groups of four or more fractures is 


P 4 - P 3 - F(CT) 


L P 3/l U - p 54/3A> +P 2/1 P 3/2 (1 - P 54/3B>’ 


(18) 


p =i - (1 .ip ) m(n - 3) 
4 8 4 


(19) 
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( 20 ) 


E 4=T P 3 M<N ' 3 > 

Although it has not been possible to obtain the general expression for the 
probability of a group containing an arbitrary number of broken elements it is 
contended that this expression is only of academic interest. It is argued that 
once the probability of secondary breaks (i. e. the probability of groups con- 
taining two or more fractures) becomes singificant failure of the composite 
can be expected. More will be said of this point later on. 


Comparison with Experiments 

The ultimate test of any theory is its agreement with experimental re- 
sults. The present theory has been compared with the data obtained in Refer- 
ences 17 and 4 for glass-epoxy composites. 

The tests reported in Reference 4 were run on composites consisting of a 
single layer of 3 1/2 mil E-glass fibers embedded in two resin systems (B and 
C). Tests were run to determine the strengths of the fibers for several gage 
lengths. As in Reference 4, it is assumed that the cumulative distribution 

function can be characterized sufficiently well by the Weibull distribution. 


F(CT) = 1 


-a Lct^ 
e 


( 21 ) 


where L is the length of the fibers and a and |3 are parameters that character- 
ize the distribution. Using this expression the variation of mean fiber strength 
with gage length is 

1_ 

= ( COL) ^ T (1+1/ g). (22) 
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The fiber parameters are found to be 


__ J_ 

P = 9 . 40 , a ^ = 181. 5 ksi. 

From the photographs of the specimens under polarized light the inef- 
fective lengths for series B and C were found to be 0. 031" and 0. 086", res- 
pectively. 

Using this data the expected number of single broken elements were 
calculated from Equation 4, and are represented by the dashed lines in 
Figures 9 and 10. The number of fractures observed experimentally in the 
various specimens in the two test series are presented for comparison. It can 
be seen that for low stress levels there are generally more fractures than the 
theory predicts. However, the behavior is most important at the higher stress 
levels in the area of failure loads, and here agreement is fairly good consider- 
ing the statistical spread in fiber properties and t he experimental uncertainties 
involved. The relatively large number of fractures at low stress levels is 
possibly a result of damage to the fibers during fabrication of the specimens 
since glass is notoriously sensitive to handling as far as strength is concerned. 
However, since the observed fractures approach the expected number of higher 
stress levels they are of little importance. If, on the other hand, failure 
occurs without a significant accumulation of fractures the breaks at low stress 
level may be of extreme importance. More will be said of this later on. 

The quantities of E , E , E and E (where E. represents the expected 
number of groups of fractures having at least i broken fibers) obtained in 
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Equations 4, 13, 17 and 20 are plotted for test series B and C in Figures H and 
12. It can be seen that the expected number of multiple fracture groups (E^, 

E^ and e 4 ) rise sharply in the observed range of composite failure. The fail- 
ure predictions of Reference 3 are presented for comparison. 

In order to assess the validity of the expressions for multiple fractures 
the number of groups of multiple fractures were counted on films of the tests. 
The results are plotted in Figures 13 and 14. The dashed curves in the figures 
are calculated values of E . It can be seen that, in general, multiple breaks 
begin to appear in the stress range predicted by the theory. Furthermore, 
the composites fail without the occurrence of a large number of multiple- 
break groups compared to the number of isolated single fractures. 

There are several factors that could account for this phenomenon. One 
possibility is illustrated by Figure 15. This graph shows the relative behavior 
of the sum P^/i + ^ 3/1 ^ or stat ^ c and dynamic overstress factors. It will be 
recalled that this sum represents the probability of the fracturing of at least 
one element adjacent to a single broken fiber. In the observed failure range 
the dynamic curve is markedly higher than the static curve indicating that 
there is a definite possibility of a failure crack being caused by the dynamic 
effects of fiber fracture rather than by cumulative static probability of 
failure. It should also be noted that the use of overstress-factors is just an 
approximation and that the actual stress concentrations caused by multiple 
breaks may be, and probably are, much more severe than those of single 
fractures. (As an analogy, the stress concentration factor for an ellipse 
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increases with the aspect ratio) Whatever the mechanism, experimental ob- 
servations seem to indicate that multiple breaks tend to occur shortly before 
composite failure. 

The present theory was also compared with two tests on continuous 

, ... , (17) 

glass fibers in epoxy run by Friedman . The Weibull parameters for the 

strength distribution of fibers used are 

J_ 

3 = 4. 0 , a ^ = 137. 1 . 

Since 3 is an inverse measure of dispersion, these fibers had a much wider 
spread in failure stress levels than did those of Reference 4 . This large disper- 
sion is reflected in the wide spread of the curves of E , E , E and E pre- 

1 2 3 4 r 

sented in Figure 16. However, failure in both specimens occurred quite 
near the stress level for which the first multiple fracture is predicted. 

On the basis of the experimental evidence cited it is proposed that the 
failure stress of a continuous fibrous composite loaded in tension parallel 
to the fibers can be reasonably well predicted by that load for which the first 
multiple fracture is expected to occur. That is, 

E 2 = WN- 1 l-( 0 ) | p 2/i+ p 3/i ) = 1. ,23) 

Analysis of Non-Cumulative Fracture Mode 
Although a large number of isolated, fractures are observed in glass- 
epoxy composites, this is not the case for other fiber-resin systems such as 
boron-epoxy and boron aluminum. For these composites failure usually 
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occurs catastrophically, without an accumulation of isolated fractures. How- 
(19) 

ever, Lenoe has demonstrated scattered fractures in a boron-aluminum 
composite containing 5% volume fraction of fibers, which is quite small. 

The absence of a significant number of isolated fractures seems to indi- 
cate that the entire composite is failing at the load at which the weakest 
fiber breaks. This would mean that the matrix is actually detrimental since 
the strength of a bundle of fibers, without a matrix, would be stronger. It 
was, therefore, decided to determine the theoretical value at which a first 
fiber fracture is expected and compare the results with experimental evidence. 

Consider a population of fibers of length L whose strength is character- 
ized by the probability density g(a). For a sample of N fibers from this group 

the distribution function for the strength of the weakest fiber has the following 
form 
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P = 7. 0, a 


_1 

3 


- 368. 


The fiber stress levels at failure for the test specimens were 193 and 215 ksi 
The expected value for the first fracture computed from Equation 26 is 171. 5 
ksi which is about 16% lower than the average failure stress of the two speci- 
mens. On the other hand, the first multiple fracture is predicted (E^ = 1) to 
occur at 300 ksi while the theory of Reference 3 predicts a failure stress of 

428 ksi. In the last two cases Friedman's definition of ineffective length was 
used. 


~ . . ( 20 ) 

urimus also ran tests on boron fibers in an epoxy matrix. These 
specimens were 2. 5 inches long and contained 25 fibers. The Weibull para- 
meters for the fibers were found to be 


3 = 11. 11, a P = 433. 


orturicitBly only one und 3 .m 3 .gecl specimen was tested. This specimen 
failed 3t 304 ksi fiber stress. The expected V3lue of stress for the first frac- 
ture is 296 ksi while the first multiple fracture is predicted at 330 ksi 3 nd the 
failure stress predicted by the theory of Reference 3 is 340 ksi. It should be 
noted that to obtain the l3st two values, the definition of the ineffective length 
presented in Reference 3 was used. 

From the experimental observation of the absence of cumulstive fractures 
3nd the good correlation between the observed failure stress levels 3nd 
those predicted for the first fiber break a good case can be made for the 
hypothesis that composites exhibiting this type of failure are only as strong 
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in tension as their weakest fiber. If this is also the case for the three dimen- 
sional composites, the matrix must be designed so that it will prevent this type 
of failure from occurring. This can possibly be accomplished by using a ductile 
matrix or by allowing the matrix to partially debond from a broken fiber. 

Conclusions 

A statistical model including the effects of stress concentrations for planar 
arrays of fibers in a matrix has been presented which provides a good descrip- 
tion of composite behavior up to the failure load. The stress level for which the 
first multiple fracture is expected to occur has been proposed as a predicted fail- 
ure stress. 

The model predicts that the composite itself is a "brittle" material in that 
its strength decreases as the length or width of the specimen increases. This is 
illustrated in Figure 17 where the variation of with length is presented for the 
composite system used in test series B of Reference 4. This is in contrast with 
the theory in that paper which predicts a composite strength that is independent 
of length for large values of M. 

Finally, it has been shown that the failure stress level in tests of three com- 
posites that did not exhibit cumulative damage occurred at about the expected value 
of stress at which the first fracture was predicted. This failure stress is lower 
than that for a bundle of fibers of the same length and number indicating that the 
matrix may have a detrimental effect in composites exhibiting this type of failure. 

The present study demonstrates that the understanding of composite be- 
havior is a continually evolving process and that even in the case of such a simple 
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loading condition as pure tension the failure mechanisms are not completely 
understood. More work must be done, for example, to explain the apparent 
difference in failure modes observed in glass and boron fibers, and to extend 
the analysis of stress concentration effects to three-dimensional composites. 
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The Steady-State Viscoelastic Response 


A systematic study of viscoelastic behavior of fiber -reinforced compo- 
sites was initiated by Hashin (Ref. 21). Such analyses are motivated by both 
theoretical and practical points of view. Fibrous composites consisting of 
linear viscoelastic phases are a natural extension of the linear elastic model 
in which phase materials are assumed to be Hookean Solids. On the other hand, 
phase materials in composites do exhibit very strong time dependent properties 
especially in a high temperature environment. A basic understanding of the 
viscoelastic behavior of composites is important before such materials are 
utilized in practical applications. 


Hashin (Ref. 21) related the effective viscoelastic properties of a compo- 
site which is composed of linear viscoelastic phases to those of its constituents 
by a correspondence principle. This same model used by Hashin is used to 
study the steady state response to some simple types of sinusoidal surface 
tractions or surface displacements. It should be noted that in Hasin's approach 
to viscoelastic theory of composites, inertia forces are neglected so that only 
quasi-static motion is treated. The present study is based on Hashin's work 
(Ref. 21) and the same assumption is implicit so that wave propagation phenome- 
na will not be considered here. 

Following Reference 21, the general macroscopic viscoelastic behavior of 
a composite can be described by the following constitutive equations: 

t 


V ‘> 3 


J G ijki 


, d e (r ) 

(t-r) kT 

dr 


(1 ) 
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or 


V‘> ■ 


t 

I 




dd, , (t) 

k-t dr 

dr 


( 2 ) 


where C7.j(t) an< ^ e -(t) are > respectively, the average stress and average strain 
in the composite at time t; G„ k ^(t) and J 7.^(0 are the effective relaxation 
moduli and effective creep compliances, respectively. 

By the symmetry of the average stress and average strain tensors, the 
following symmetry relations hold: 


G Jikl (t) ' G ’ijM <t) = >..«> 


for t > 0 


(3) 


J jik^ (t) ' J ijk^ (t) J ijtk (t) 

In order to establish the correspondence principle between viscoelasticity 
and elasticity, it is assumed in Reference 21 




and 


for t >0 


(4) 




which are the Onsager Reciprocal Relations (Ref. 22) 

The one-sided Laplace transform of (1) and (2) gives 


V p) = p G uk-c <P> V<P> 


(5) 
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( 6 ) 


V p) = p J !jk<. (p) v (p > 


where p is the transform variable and the circumflex " / " above a function 

denotes its Laplace Transform which is assumed to exist. Because of the 
formal resemblance of (5) and (6) to the generalized Hooke's Law in elasticity, 


p and PJ'." Hp) 

ijk' 0 ijk'l are terme 

compliances, re spectively. 

Now for unidirectional fiber 


d transform 


-reinforced 


domain (TD) effective moduli and 


composites, transverse isotropy 


will be considered (Ref. 21 , 11 ). Consequently, only five effective relaxation moduli 
(or five effective creep compliances) are independent so that (5) and (6) can be 
much simplified. For example, (5) can be written down in terms of the following 
five independent TD effective relaxation moduli: 

P ^^(p) plane strain TD effective bulk modulus; 

a , 


P g 23 (p) - 
P Gj (p) -- 
P E* (p) -- 


the TD effective transverse shear modulus; 
the TD effective in-plane shear modulus; 
the TD effective Young's modulus; 


Aj|; 


? C n(P)-" to be associated with uniaxial stress in fiber direction 

with transverse deformation prevented by a rigid 
enclosure. 

On the other hand, the TD effective creep compliances are the reciprocals 
of the corresponding TD effective relaxation moduli. For example 
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1 


P k 2 3 (P) = a * 

PK 23 ( P ) 


/ 

P g 


23 


(P) 


p g 2 * (p) 


etc. . 

After the effective characteristic functions (relaxation moduli and creep 
compliances) are defined, Hashin(R ef.2l) used a correspondence principle to relate 
the effective viscoelastic characteristic functions to the effective elastic moduli 
of a duplicate composite body with elastic phases. 

On the other hand, the effective elastic moduli of unidirectional fiber- 
reinforced composites have been obtained by Hashin and Rosen(R e f. 11). Closed-form 
expressions for four effective elastic moduli and bounds for the fifth are ob- 
tained for the "random array" geometry. In applying the corresponding prin- 
ciple, Hashin (Ref. 21) used these expressions and went through a replacement 
scheme to obtain the Laplace Transforms of the corresponding effective visco- 
elastic characteristic functions. Therefore, the entire problem of finding the 
quasi-static viscoelastic response reduces to Laplace Transform inversion 
which is not always easy. However, for cases where inversion is formidable, 
Abel-Tauber theorems can be used to draw important conclusions on the be- 
havior of the effective moduli and compliances without the operation of inversion. 
Refs. 21 andll are referred for details. 

Now, for a unidirectional fiber-reinforced composite body under special 
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where a (t) ( e (t) 

whereas G := (t) represents one of the five effective relaxation moduli 

(effective creep compliances) of physical importance as defined previously. 
Henceforth, (7) (or (8) ) will be used as the representative average str e s s - strain 

relation that defines a particular viscoelastic characteristic function. 

In the course of obtaining the steady- state response by making use specifi- 
cally of the results of Ref. Zl, we shall first develop certain relationships and 
investigate their validity. Then the general results will be applied specifically 
to the unidirectional fiber- reinforced composite to obtain explicit expressions 
which characterize the response. Since we are interested in steady-state re- 
sponse to sinusoidal input, we can assume that the input has been applied on the 
body for an indefinitely long time and that all initial disturbances have died out. 
Under this circumstance, it is convenient to put the beginning of motion at 
time - 00 (Ref. 2 3). Hence (7) and (8) are modified to the following forms: 

5(t)= j G*(t-T) d - e J T - QT (9) 


is a component of the stress (strain) tensor 0\. (t) (e..(t) 
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and 


e (t) = J J (t-T) 


d g (r) dr 
dr 


( 10 ) 


By changing the integration variable from T to £ where t - T = £ , ( 9 ) and 
(10) become, respectively, 


a(t) = t G (t) V (t-i)di 


/ 1 1 \ 


1 


e (t) = i ra)o (t-4)d4 

o 


(12) 


where prime denotes differentiation with respect to the argument of the function. 
Using complex representation for sinusoidal oscillation, we put 

iCOt 


e (t) = e e 
o 


( 13 ) 


where is in general a complex number and 60, a real number, is the angular 
frequency of the oscillation D 

Substituting (13) into (11) we obtain 


- i 60 1 

O (t) = icoe e G (co) 
o 


( 14 ) 


where 


G*(co) j G*U) e" lU ^d£ 


( 15 ) 
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is the one-sided Fourier Transform of G (t) if it exists. 


itOt 


Let Cf(t) - CT^e , where or is in general a complex number, then by (14), 


CT(t) a o • * 

= -2- = 1 OJ G (co) M 
fi (t) C 


( 16 ) 


where M'“ is called the effective complex modulus associated with G(t). 
Equation (16) can also be written as 


cr(t) = Im'I e e 
1 o 


i (cot + 6 ) 


(17) 


where 6 = tan 


-1 


Im [M ] 
Re [m"] 


which reveals the fact that O (t) is not in phase with e (t) 


Similarly, if a (t) = a^e lC ° t then (12) becomes 


€ (t) = i COCT e 1 ^ J ' (CO) 
o 


(18) 


where 


J*(U) = j J* (i) e -U °^d i 

o 


(19) 


is the one-sided Fourier Transform of J '" (t) if it exists. 


Let e(t) = e e 
o 


icot 


then by (18), 


l (t) _ w j* (co) = 

5 (t) ° 0 M" 


( 20 ) 
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or 


Mt) = ^C e i(Wt - 6 > 
M o 


( 21 ) 


which again shows that e (t) lags behind CT (t) with a phase lag & . 

From both (16) and (20), it follows that 

2 ''"'ye 

- CO G"( CO ) j !,: (co) = l 
if both G (CO) and J (co) exist. 

The physical significance of Im [ M ] and R e [M ] is now discussed. 

First of all, it can be shown that the rate of work done by surface traction at time t 
on a body can be expressed in terms of an integral over the entire body as follows: 

I T. u. d S = I a 

•iii J n 


e ,, d V 


ij ij 


where S denotes the boundary surface and V, the volume of the body. The above 

equality is valid only under the assumption of quasi- static motion. 

Furthermore, under special boundary conditions on the boundary surface 

S (displacement boundary condition u.(S) = e..(t) x. or traction boundary condition 

i ij J 

T.(S) = cr..(t) n ,(Ref. 21)), it is easy to show that: 

1 !J j 

l cr. . (x, t) e.. (x , t) d V = a..(t) e.. (t) V 

J iJ “ b ~ ij ij 

V 


Therefore, the total work done from time t^ to time is: 


J L ) T i dS 


f _ _L 

dt = V \ a..(t) e .. (t) dt 

J iJ ij 


t 


1 
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If t - tj is the period of motion, then 


W = V 


f 


2 5 . (t > €..(t) dt 

1J IJ 


( 22 ) 


denotes the total work done on the system in a cycle* According to the First 
Law of Thermodynamics, for isothermal steady state deformation, W is the 
total energy dissipated and transferred to the surroundings in the form of heat 
in a cycle. Under the special boundary conditions by which (7) and (8) are ob- 
tained, W in (22) can be reduced to the following simple form: 

L 


w = VF ( a (t) e (t) dt 

J 

h 


(23) 


where the factor F is either 1 or 2 depending on the boundary condition. 

iCOt 


Therefore, if e (t) = Re [ e e ± '* Jl ' ] and CT(t) = R [ivT'e e lCOt 

o e o 


] , then for 

a cycle with angular frequency of motion CO, after some manipulation, (23) 


becomes ZlJ 

T + 


W = VF \ 5 e dt = 77 e Im [ M ] VF (e is assumed 

o o 

T 


1 


to be real without loss of generality ) where T is any time during the motion. 
This gives the amount of energy dissipated in a cycle. 

Furthermore, it can be shown that a fraction of the amount of energy 


° |R e L M'j I Vb is twice in every cycle alternately stored and expended in 
2 
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the system. Therefore, Im [M‘ ] is related to the energy dissipated while 

Re [ M r ]is related to the energy stored. 

From the above results, it is clear that in the steady state, the 

responding average stress (or average strain) will vary sinusoidally with the 

same angular frequency CO as that of the input average strain (or average stress) 

with a phase difference & which is a function of CO . Moreover, the amplitude 
a 

ratio is also a function of co only. 

G o 

According to (16) and (20), the knowledge of M hinges on the knowledge 

' * ' 5*C 

of G (co) or J (co) . However, G (co) and J (co) are formally related 
to the Laplace Transforms G"(p) and J"(p) by the following relations: 


and 


G* (CO) = G*(ico) 


J v ( co ) = j'(i to) 


(24) 

(25) 


Therefore, if expressions for G (p) and J (p) are known, G (co ) and J (co) 
can be obtained immediately by the replacement of variable p in G (p) and 
J (p)byiUj. However, there still exists the problem of convergence -- 
given a function of time t, the existence of its Laplace Transform does not im- 
ply the existence of its Fourier Transform. Thus, given a Laplace Transform 
of a viscoelastic characteristic function, it is necessary to examine the location 
of its singularities in a complex plane of complex variable p. According to the 
definitions of the one-sided Laplace and Fourier Transforms, it is observed 
that if the singularities of the Laplace Transform are all located in a region to 
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the left of imaginary axis of the complex plane, then the corresponding Fourier 
Transform exists. 

Hashin (Ref. 21) has obtained closed-form expressions for some effective 

relaxation moduli and creep compliances for the "random array" model. Then, 

by making use of (24) and (25), the corresponding effective complex moduli can 

be obtained in a straight-forward manner. For example, for elastic fibers 

embedded in a viscoelastic matrix, the effective relaxation moduli K' (t) 

2 3 

characterizing the plane strain dilatation has the following form in the Laplace 
Transform domain: 


P k !(p) = [k ( P ) +~r ( P )] + v 

£ 3 m 3 m f 


m 


- 1 


K„ 


+ 3 C f- 


k J P > 


3 r J»> 


k m <P)+ I r m <P » 


where k^p) , ^(p) are matrix TD moduli; K f and G f are fiber elastic bulk 
and shear moduli; v^ and v^ are the fiber and matrix volume fractions, respec- 
tively. All these quantities are considered known if the composite body is 
given. Then according to (24), formally we have: 


m 


icuK' (w) = [k (ico) + ^ T ( ico) ] + v 
^ J m 3 m f 


K f + 3 G f 


k (ico) - - T (ico) 
m 3 m 


k ( ic0 ) + T r (i to) 

m 3 m 


and the complex moduli associated with K (t) is, according to (16) , 

Lt Jj 


M 23 = iC ° K 23 (W) 


Furthermore, it follows from Ref. 21 that for rigid fibers. 


K 23 ( ")= K m <“» + 5m<“) 

3 


’+ [^(0.)+ | G m < “ ) ] (1-v) 
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where K^(co) and G^tO ) are the one-sided Fourier Transforms of the matrix 
bulk and shear relaxation moduli, respectively. On the other hand, if fibers 
are rigid and the matrix is elastic in dilatation and Maxwellian in shear, 
Hashin(R e f. 21)has obtained simple expression for the effective creep compliance 
k 23^) which the one-sided Fourier Transform does not exist. Similarly, 
other steady- state responses such as in- plane shear, etc* can be obtained in a 
straight-forward manner. 


It is emphasized here that the method developed here is also valid for 
other models than the composite-cylinder-assemblage model used in Ref. 21 ; th 
latter is only a case in which bounds of some of the characteristic functions in 
transform domain coincide to yield closed-form expressions. If other ex- 
pressions can □ e o o Laine u uy other models, the corresponding effective com- 
plex moduli can be obtained by the replacement scheme defined by ( 24 ) and (2 5), 
On the other hand, vibration experiments can be performed on composites to 
determine experimentally their respective effective complex moduli from 
which information on their corresponding effective viscoelastic characteristic 
functions can be obtained through Fourier transformations. 
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Evaluation of Transverse Effectiveness Factors for Use in Elastic Analysis 
of Three Dimensional Filamentary Composites 

The method of analysis developed in Reference 1 for the elastic constants 
of composites having filamentary reinforcements at various angles to the three 
principal orthogonal axes of the material employed factors f ) P • * 


etc. to define the transverse effectiveness of the filaments for resisting stretching, 
sneanng, etc. For simplicity the assumption was made that the transverse 
effectivenesses were the same for all strains due to stretching ~ regardless of 
whether the strains were induced by Poisson's ratio effects, or were simply 
those in the direction of an applied extension. (Different values of transverse 
effectivenesses in shear from those in stretching were, however, allowed for.) 
inis approximation led to generally satisfactory results for stiffnesses, with 
the possible exception of E^, that is the stiffness along the filaments, for uni- 
directional reinforcement - a somewhat disturbing exception in view of the 
simplicity of calculation of this stiffness by the rule of mixtures. Even more 
disturbing was the fact that the values of Poisson's ratios calculated with this 
simplifying assumption of equal transverse effectiveness factors were not con- 
sistent with those found by other methods of calculation. 

line s for the extension of the analysis of Reference 1 to provide for 
various values of ^ for direct and Poisson strains were included with the 
analysis, but within the time available no evaluation could be made to determine 
whether or not the extension would be justified. A brief series of calculations 
has now been made to evaluate the differences among the various transverse 
effectivenesses for the extended analysis. The results are shown in Figure 18, 
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and the equations relating thereto are reproduced as Tables 1 to 3. 

In Figure 19, curves of 9 G are aS 

derived from the upper limit values of the elastic constants of Reference 11 for 
a typical glass/ epoxy combination employing the assumption that 

ft ^ can seen, differences among 

, and are found, of increasing relative magnitude with in- 

creasing volume fraction of binder (the abscissa on the Figure). Also 
differs from the value unity; this is the variation which is to be anticipated 

from the previously noted inconsistent calculation of E for uni- dir ectional 

-Lj 

reinforcement with the simplified transverse effectiveness assumption. The 
differences in the various - values shown are, of course, just those consistent 


with the elastic constants as found from F.efererme 


The use of the A of 

r ’ 


Tables 2 and 3 thus reproduce properly all the elastic constants 
of a unidirectional reinforcement configuration, and hence provide a 
self-consistent basis from which the effects of multi-angular reinforcement 
may be determined via extensions to equations like those of Reference 1, 

These extensions have been made and the general equations for the compliances 
incorporating the various are presented in Table 4. 

As presently written, the equations of Tables 2, 3, and 4 conceal some 
internal inter-relationships among the various expressions. For example, the 
trigonometric expressions of Table 4 are for convenience written in terms of 
the three angles ^ , fjf , and which the reinforcing filaments make 

with the composite orthogonal axes of symmetry (1, 2, and 3). Only two 
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angles are needed, however, to define the filamentary directions; 
and ii are related by the well known formula 



COS*f + c Ot l j> * COS 1 A s / (1) 

Similarly, for convenience six transverse effectivenesses have been 
defined } /$£ } ft', ^ see tables 2 » and 3 )> 

where only five are needed for consistency with the five elastic constants 

used in evaluating the transverse effectivenesses. Accordingly the equations 
of Table 3 may be combined, as for example to yield the relation 



'V 


(2) 


where u = the Poisson's ratio of the filamentary material. 

Thus simplifications of the equations as given in Table 4 are undoubtedly 
possible by the employment of the trigonometric expressions, and values 

which lead to the least complex algebra. 
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DESIGN CONCEPTS 


Study of Composite Structu res 

Efficient application of composite materials to aerospace structures 
requires proper selection of reinforcement pattern and material as well as 
overall structural arrangement. To a large degree guidelines for optimum 
design of such composite structures have been lacking, nor have they been ob- 
vious a priori to the designer. 

Some clues about possible directions toward efficient configuration have 
recently become available. Reference 24 explored effects of material and rein- 
forcement pattern on the structural efficiency of boost vehicle shells. Reference 
2 investigated similar effects for a number of aircraft structural elements. 

While these studies revealed a number of specific factors of importance for 
efficiency of application of composites, three general conclusions also evolved 
which served as a basis for the studies to be reported here. These conclusions 
were: 

1. For a wide range of shell type applications the isotropic (0°, ±60°) 
reinforcement configuration is most efficient. 

2. The high (multi-directional) stiffness-to-density ratio of beryllium 
makes it more efficient than most near-future composites for many 
shell and plate type structures, - particularly if the load intensities 
encountered are low, or if stiffness requirements are important. 

3. Filamentary composites appear most attractive when used as 
unidirectionaliy reinforced elements to carry unidirectional loads. 
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The further question that seems to arise from these conclusions is 
whether some combination of the biaxial properties of beryllium (or isotropic 
composite) with the uniaxial properties of one-directional filamentary rein- 
forcement may possibly be the best configuration of all. 

Circular Tube -Columns 

To explore this possibility on an orderly basis, a simple round-tube 
column was selected as a first model for study. This model was chosen for 
simplicity; for an isotropic material optimization of the thin-walled tube in 
compression has perhaps achieved the status of being a classic example of 
balancing proportions between local buckling and column bending to achieve 
minimum weight. Thus its use to explore effects of combinations of unidirectional 
and multidirectional properties is uncomplicated by complexities of the optimiza- 
tion procedure itself. 

Accordingly beginning with the model of Reference 25 , modifications were 
assumed of increasing complexity, and the effects on efficiency evaluated. As 
a first step, unidirectionally reinforced composite stiffeners were assumed added 
to the tube at three equally spaced points around the circumference. These 
stiffeners were assumed to contribute nothing to the local stability,-or instability, - 
of the tube walls, but were assumed to add effective areas for column bending 
so that the effective radius of gyration of the tube is given by the expression 
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where 



Eff 


R 


t 


E 


B 


E 

o 


v 


B 


v 

o 


radius of gyration of stiffened tube 
outside radius of tube 
thickness of tube 

Young's modulus of unidirectional (boron) composite 
loung's modulus of tube material 
volume fraction of composite stiffeners 
volume fraction of tube material ( v + v =1) 


With equation (1) incorporated into the procedure of Reference 25, the 
efficiencies of the reinforced tubes were calculated from the formulas that 


•foil 


P_ 


- 0.00912 



( 2 ) 


K. fryfesr 
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where 


axial load 
column length 
overall weight 

density of stiffening material 
density of tube material 

and the subscripts Sec and Tan refer to the secant and tangent moduli of the 
tube material, respectively. 

As indicated in the formulas, the reduced moduli used for plastic buckling 
were the tangent modulus for column bending, and the secant modulus for local 
buckling. The use of the tangent modulus for column bending is well founded. 
The use of the secant modulus for local buckling is used here as slightly more 
optimistic than the root mean square of the tangent and secant moduli as pro- 
posed in Reference 26. Thus the efficiency curves resulting are perhaps slightly 
too high in the plastic region, representing a kind of upper bound. (Most 
affected are the curves for beryllium which should be accordingly somewhat 
discounted at the upper end.) In any event the use of this possibly optimistic 
reduced buckling modulus changes none of the conclusions drawn from the 
results. 

Results of the calculations are given in Figures 1 9 to 24. In Figure 19 are 
given basic results for unstiffened aluminum-alloy and beryllium circular-tube 
columns. Here the beryllium tubing is found to vary from 2 5% of the weight of 

P/ L* K % C 0 -*'*'- 


the aluminum-alloy tubing in the elastic range 


■ 



p 

L 

W 
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to 80% of the aluminum-alloy weight at high stresses (%>** *£ [*op*n) 
The addition of the boron/epoxy three-point reinforcement permits the 
aluminum-alloy tubing to be made substantially lighter both at high and low 
loadings (Fig. 20). In the elastic stress range there is apparently an optimum 


reinforcement ratio 


(£* °- 667 ) 


beyond which higher percentages of 


reinforcement do not further increase the efficiency. At the high stresses, 
however, the very great strength and stiffness of the unidirectional reinforced 
composite provide increases in efficiency up to the maximum reinforcement 
ratio considered 



The beryllium round tubing is not improved by the boron composite rein- 
forcement at low stresses; rather the efficiency is decreased as the reinforce- 
ment ratio increase S: This trend is just barely reversed, as might be expected, 

at high loadings (see Fig. 21). 

Comparisons between Figures 20 and 21 reveal that at low loadings the 
unreinforced beryllium is always the lightest of the combinations considered, 
being approximately one -third the weight of the best boron/epoxy- aluminum. 

At the high stresses on the other hand the reinforced aluminum is better than 
any beryllium tube, being about one -half the weight of the beryllium construc- 
tion at the maximum loading considered. Thus it appears that the low propor- 
tional limit strain of the beryllium together with the flat top to the beryllium 
stress -strain curve beyond the proportional limit prevents it from acting 
effectively in conjunction with high-strength uni -directional reinforcement. 

In this range a better combination of materials is an isotropic boron/epoxy 
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configuration for the tube walls. This configuration should still be elastic at 
the values of P / ^2 covered in Figures 19 - 21, and at the maximum values of 

-p 

/ l 2 considered the isotropic boron/epoxy tube is calculated from equations 
(1) and (2) to weigh very nearly one -half as much as the best reinforced 
aluminum tube. 
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Scalloped Tube Columns 


One of the pitfalls of efficiency studies like the foregoing is the possibility 
that a poor geometrical shape has been chosen, and as a result misleading con- 
clusions are derived about the effects of various materials of construction. 
Ideally an optimum shape should be used (if need be the optimum for each 
material combination) so that shape effects may be divorced from material 
effects in the evaluations of the calculations. 

In order to assess the shape effect, the efficiencies of scalloped tubes 
were also investigated. First, a three -lobed scalloped round shape was used 
(Figure 22). The objectives of the scalloped shape are enhanced local buckling 
resistance due to decreased radius of curvature of the tube wall together with 
increased column bending resistance because the. scalloped tube approximates 
an equilateral triangle, and the equilateral triangle has a 21% greater cross- 
sectional moment of inertia for the same area than a circle. Because the 
scalloped round tube does not substantially enhance the moment of inertia, 
however, it was abandoned in favor of a scalloped triangular tube. The charac- 
teristics of these two shapes are summarized in Figure 22. 

In Figure 22 are plotted the percentage increases from a simple thin-walled 
round tube section attained by scalloped-round and scalloped-triangle sections 
having the same cross-sectional area and radius of curvature -to -thickness 
ratios as the reference round tube. As indicated by the curves of the Figure, 
the "scalloped triangle" has substantially greater potential for enhancing the 

T> 

column strength (at a constant / t so that in first approximation the local 
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buckling strength is constant) than the "scalloped circle." Furthermore, the 
apexes of the triangular shape provide a greater radius of gyration for added 
stiffening as in the form of concentrated boron/ epoxy elements than the crests 


of the scallops of the "scalloped circle." 

While clearly the extreme limit of 
represented by the "infinitely scalloped" 


21% increase in section properties 
triangle f&sOJ of Figure 4 surely 


(f-0) 


cannot be attained, the development of a reasonable fraction of that amount may 
be anticipated for "reasonable" proportions. Here "reasonable" is taken to 
mean that the cusps between scallops are of sufficient depth to establish stable 
corners not prone to local buckling at the stresses encountered by the tubes. 

The assumption was made that the proportions so labelled on Figure 22 are 
reasonable, and likely performance gains through their use for reinforced com- 
posite columns was calculated with the results given in Figures 23 and 24. For 
these calculations, equation ( 2) was replaced by 



/h« 


re 



moment of inertia of scalloped triangle section 
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area of triangle section (equal to area of circular 

■p 

section of equivalent / t ) 
inside radius of scallops 
equivalent thicknes s -radius ratio of tube 

In Figures 23 and 24 comparisons of the efficiencies of the reinforced 
scalloped triangle tubes with envelope curves representing the lightest reinforced 
round tubes (the dotted curves on Figures 23 and 24, derived from Figures 20 and 21) 
show the expected gains for the triangular shape. With the triangular shape 
higher reinforcement ratios for the aluminum alloy tubing are effective (Fig. 23) 
and a reinforcement ratio = 4 produces a tube column competitive with 

beryllium even in the elastic range (see Fig. 21). Further, the boron/epoxy re- 
inforcement of beryllium in the plastic range is more effective with the scalloped- 
triangle shape. The general trends established in the studies of the remforced- 
round tube columns, however, are not changed by the more efficient triangle 
shape. Thus the scalloped-triangle beryllium tube is the most efficient of all at 
low loadings and the reinforced aluminum -alloy becomes more efficient when the 
stresses would cause yielding of the beryllium. 

The gains in efficiency shown on Figures 23 and 24 for the scalloped-triangle 
shape appear great enough so that experiments to determine their validity appear 
desirable. That is, experimental definition is needed of the stability of a 
scalloped-wall tube to establish the degree of scalloping beyond which corners 
of cusps will no longer remain fixed as straight lines along the tube to provide 
edges as restraints against local buckling distortions. 
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Plate Efficiencies 


The efficiencies of composite materials for plate applications, already 
investigated for a variety of material combinations in Reference 25, have 
exhibited similar characteristics to those already discussed above for tube- 
columns. Indeed combinations of beryllium and 0° boron/epoxy plates made 
up as sandwiches with the beryllium as the faces and the 0° boron/epoxy as 
the core demonstrate rather clearly the problems encountered in the develop- 
ment of combination metal composite structures, as is shown in Figure 25. 

In Figure are plotted the results of calculations (made as in Reference 
25) of the efficiencies of sandwich plates which have 0° reinforced boron/epoxy 
cores and faces of beryllium, isotropic boron/epoxy, and (at the very top end of 
the curve) no faces at all. The solid curve represents the beryllium -fac ed 

AW 
m* 

are totally elastic (to achieve this elasticity at such high values of a com- 

pressive pre -stress is assumed in the boron/epoxy core just sufficient to stress 
the faces to the proportional limit in tension at zero external load). The rapid 
increase in weight of sandwich above the discontinuity arises from the continuing 
decrease in stiffness of the beryllium faces as they are stressed further and 
further beyond their proportional limit. The dotted curve represents the higher 
efficiency that would be attained with beryllium with a higher elastic limit (or 
pre -stress). Inasmuch as neither of these improvements in the beryllium 
response appear possible of attainment in practice, the solid curve is the 
realistic one for the sandwich, and, as indicated, it can be surpassed in 


[iJ Msi] ) the plates 


N x 

sandwiches; up to the discontinuity (at 
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efficiency - as by sandwiches utilizing isotropic boron/epoxy faces at the higher 
load intensities. Indeed as the loads become high enough, only the 0° boron/ 
epoxy retains enough stiffness to be effective, and no faces at all become most 


efficient. 


To a degree, the curves of Figure 25 may be considered to depict in 
general the various characteristics which must be balanced in composite plate 
construction for maximum efficiency, as follows: 

(1) 0° reinforced material should be used as the strength element, 
if possible pre -stressed to provide the maximum possible 
elastic range to the material used to provide transverse 
stiffness or plate buckling resistance. 

(2) The material used to provide the plate buckling resistance 
(the face material in a sandwich) should be selected on the 
basis of the stresses to be carried, — at low stresses beryl- 
lium is most attractive, at higher stresses composite con- 
figurations approaching closer and closer to the 0° config- 
uration should be used. 


(3) While the example shown (Figure 25)is for a solid sandwich, 


so that the average stress is simply 


A£_ H,/b) 

t 'bit/ 


, the 


same rules as (1) and (2) apply for hollowed out constructions 


which raise the average stress to ^x where t is the aver- 

T 

age thickness and T < t . Thus, in general, any hollowing 


out (light weight core sandwich construction, etc.) which 
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does not introduce new instabilities or failure modes is 
like an increase in the structural index value on Figure 25 f 
- i.e. , it tends to make optimum constructions utilizing 
higher percentages of 0° reinforcement and lower 
percentages of beryllium. 
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Panel Efficiencies 


The plate efficiencies just discussed were considered on the basis of 
initial buckling stresses only. Particularly at low stresses, as is well known, 
plates of most materials can continue to carry compression load without failing 
at stresses substantially above the initial buckling stress. 

For plate assemblies incorporating 0 reinforced composites as stiffeners, 

initial buckling of the plates may be expected to cause the major part of any 

subsequent load increase to be borne by these stiffeners. This load transfer 
o 

into the 0 reinforcement is precisely the mechanism noted desirable above for 
increased structural efficiency. Accordingly, for example, compression 
panels incorporating 0° reinforced stiffeners should be expected to achieve 
high structural efficiencies. 

High structural efficiencies were indeed found for 0° boron/ epoxy Z- section 
stiffeners on ±45 reinforced boron/ epoxy skin in Reference 2. Part of this 
high efficiency arose from the high shear stiffness of the skin material 
(allowing the use of thin skins and hence allowing most of the material to be 
incorporated in the stiffeners.) 

Reviews of the properties of ±45 boron/epoxy and comparisons with the 
properties of beryllium suggest that replacement of the ±45° boron/epoxy skin 
on the panels of Reference 2 with beryllium should result in even higher effi- 
ciencies than those calculated therein. For example, pertinent properties of 
these two materials are compared in the following table: 
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±45 B/Epoxy 

Be 

Density 

2. 19 

1.86 

(pci) 

(0. 079) 

(0.067) 

MN 

Yield Stress — — 
rm 

138 

400 

(ksi) 

(20) 

(58) 

Shear Modulus 

72.5 

145 

m^ 

(ksi) 

(10, 500) 

(21,000) 

Young's Modulus GN 

25. 5 

304 

2 

m 

(3700) 

(44, 000) 

(ksi) 



Yield Strain 

0. 0054 

0. 0013 


Thus, except for its low yield strain, beryllium appears vastly superior to the 
i 45 B/Epoxy. The low yield strain, however, raises uncertainties, and a 
detailed evaluation is required to determine whether the beryllium/ composite 
panels are as superior as would at first appear. 

Such an evaluation was made as a part of the study of Reference 27. The 
results are reproduced here and extended to make possible direct comparisons 
of the efficiencies of the ±45 B/Epoxy and beryllium- skin Z-panel constructions. 
These comparisons are presented in Figure 26. 

Figure 26 plots curves of weights of box-beam compression covers (plus 
supporting ribs) for optimized Z-stiffened panels. Optimization includes the 
selection of rib spacing and stiffener size and spacing for minimum weight to 
carry a design bending moment M with a box beam of width b and depth d. 

Further the skin thickness is required to be adequate, in terms of the box depth, 
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to provide a specified torsional stiffness; - accordingly the beryllium and 
±45 B/Epoxy are compared at "equivalent" values of skin thickness t 


such that 


Q/ipoxy 




B/epoxy 


)- 


2 t s, 


The curves of Figure 26 show substantial weight savings for the beryllium 
o 

skin on 0 B/Epoxy stiffeners. These weight savings are depicted by the cross- 
hatching between the curves for the ±45° B/Epoxy-skin panels (the dashed 
curves) and the beryllium- skin panels (the solid curves) at the equivalent tor- 


sional stiffness measures (equal values of 


d A 


Several factors accumulate to produce the weight savings shown in Figure 
26 for the beryllium- skin composite construction: (1) the favorable shear stiffness 


of the beryllium permits an even higher percent of stringer mater ial-with its 
high-strength, high column-bending stiffness characteristics, - than for the 
±45 Boron/Epoxy skin; (2) the beryllium is of lower density than the boron/ 
epoxy; (3) the high strength of the 0° reinforcement (as noted in the opening dis- 
cussion for this section) permits it to carry the load long after initial buckling 
of the beryllium skin. 
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Studies of Composite Materials 

Three -Dimensionally Isotropic Materials 

The form of the matrix of elastic constants relating stress to strain in a 
generalized Hooke’s law has been studied for various conditions of structural 
symmetry. For two dimensions, it is known that a material is isotropic in a 
plane for which there exists a normal axis of at least six-fold symmetry. Thus, 


a s 


a practical example, a fibrous composite plate is effectively isotropic in 


its plane when one n th of the fibers are oriented every -1 radians for n » 3 

n 

(By effectively isotropic, it is meant that the average stress, average strain 
relations are isotropic. ). Similar conditions to obtain three dimensional iso- 
tropy appear to be unavailable in the literature. Several configurations having 
multiple symmetry conditions which might yield at least an approximation to 
an elastically isotropic material have been studied. The first material treated 
has two three-fold axes of symmetry separated by the angle for which the cosine 
is -1/3. Such a material can be obtained with a fibrous composite having one 

fourth of the fibers oriented in each of the four directions defined by the altitudes 
of a regular tetrahedron. 

The stress-strain relations for a general elastic body may be written as: 


T. = c. e 

1 D j 


i. j = 1, 2. . . 6 


(1 ) 


where j are the stress components 

G j are straan components 


C ii are t * ie ^ independent elastic constants (c = c ) 

J ji ij 
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We treat a material having two axes of three-fold symmetry at an angle 0 


such that 

cos 0 — — 

As a consequence of this, there must exist two additional axes of three-fold 
symmetry; the four axes being oriented along the directions of the four alti- 
tudes of a regular tetrahedron. (This four-legged array is called a caltrop. ) 
We select the vertices of the tetrahedron at the points: (0. 0. — ), 

{ W ' °’ "tJ~ T ) ’ 


( - 1 — J_ [Z) t 1 1 _L IT) 

2/3 ’ 2 ’ 4 ** 3 ’ ' 2/3 2 ’ 4 V 3 ‘ 


This is a tetrahedron with centroid at the origin. This body has elastic symme- 
try with respect to a rotation of — ^ about the z axis. Thus for the following 
transformation of coordinates, the stres s - strain relations remain unchanged: 



X 

y 

z 

x' 

COS 0 

sin 0 

0 

y' 

- sin0 

cos 0 

0 

z' 

0 

0 

1 


This reduces the elastic constant matrix to the following form (Ref. 28 ) 


ij 


c 

11 

C 12 

C 13 

C 14 

c 

15 

0 

C 12 

c u 

C 13 

- c 

14 

- C 15 

0 

C 13 

C 13 

C 33 

0 

0 

0 

i. j. 

" C 14 

0 

C 44 

0 

_C 15 

O 

1 — * 
(J1 

’ C 15 

0 

0 

C 44 

C 14 

0 

0 

0 

• C 15 

C 14 



( 2 ) 
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This array has seven independent constants. 

The stress-strain relations will also be invarient for a rotation about 
one of the other axes of three-fold rotational symmetry. In particular, we 
consider a 120 ° (counter-clockwise) rotation about the axis passing through 
the point | \ followed by a 12 0 ° (clockwise) rotation about 

v // 1 ,rz\ 

the axis passing through the point/ -p~- , 0, - — J — j . The direction cosines 


for this transformation are 



X 

y 

Z 

x' 

1 

3 

0 

2/2 

3 

y' 

0 

-l 

0 

Z 1 

2/2 

o , 

1 


It can be shown that the elastic constant matrix must therefore be of the 


form 


C 12 

C 11 +C 12' C 33 

0 

/ 2 (c ir c 33> 

0 

C 11 

C 11 +C 12 _C 33 

0 

-/2( c ii-c 33 ) 

0 


C 33 

0 

0 

0 



3 

2 C 11 " 

C i 2 

— -c 0 

2 33 

- /2(c ir c 33 ) 


T C 11 2 "~ C 33 ° 

_ T ,C u' C 12 ) 

This material has only three independent elastic constants and is a 
material with the symmetry of the cubic crystal. Indeed, a transformation 
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of the principal axes of a cubic crystal given by: 



X 11 

y" 

z" 

x' 

l 

n 

l 

72 

0 

y' 

i 

7 6 

i 

76 

-J7 

Z 1 

1 

n 

1 

77 

i 
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followed by a transformation of axes given by 



x' 

y' 

z’ 

X 

0 

-l 

0 

y 

1 

0 

0 

Z 

0 

0 

1 


brings the principal axes of a cubic crystal in coincidence with the x y z axes 
which are the reference axes for (5). This transforms the array of elas- 

tic moduli of a cubic crystal to the following form; 


c.. 

ij 


^-(C+D+2E) 


(C+5D-2E) 

6 

y(C+D+2E) 


y(C+2D-2E) 
-i- (C+2 D-2E) 
Y (C+2D+4E) 


0 T7? (c_d “ 2e) ° 

0 — — (C-D-2E) 0 

0 0 0 

f( C - D + E) 0 ^(C 

y(C-D+E) 0 



-D-2E) 

D+4E) 

(4 ) 
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where C, D and E are the three independent elastic moduli for the cubic 
system. 

The array (4) has the same form as the array (3). In fact, if we let 


C = 4c ll- 3c 33 


D = 


•c. + 6c + 6c 
“ 12 33 


and 


E = 


12 

(C U +C 12’ 


•4 c 


33 


then (3) and (4) are identical. This shows that the material with the elastic 
moduli of the form of (3) is elastically cubic. 

Although this material is not isotropic, the properties for a particular 
fibrous composite prove interesting, as will be shown in the following sec- 
tion. 

Another configuration of interest is that of Reference 1 which is construc- 
ted by taking three pairs of reinforcing filaments, oriented with reference to 
an orthogonal Cartesian coordinate system, xyz, as follows: one pair in 

the xy plane making angles of + 0 with the x axis; one pair in the yz plane 
making angles of +0 with the y axis; and one pair in the zx plane making 
angles of j-0 with the z axis. It was shown in Reference 1 that when 0 is ap- 
proximately 30 , the shear modulus, Poisson’s ratio, and Young's modulus 
satisfy the relation: 

c - E 

2 ( 1 + 1 /) 
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It can be shown that when 


0 = tan (2 sin 18°) ~ 31° 43' 

the twelve rays ( 6 lines) from the origin pass through the vertices of a 
regular icosahedron with centroid at the origin. These six lines are there- 
fore axes of five-fold symmetry. The possibility that a body having six 
axes of five-fold symmetry as isotropic is suggested. This result is per- 
haps of academic interest only, as it does not appear possible to construct 
a continuous space lattice having five-fold axes of symmetry. 


85 



Evaluation of Moduli of Caltrop Reinforced Materials 
The equations derived for the elastic constants of three -dimensionally- 
reinforced composites were used to compute properties of the caltropic re- 
inforcement configuration. Calculations were made using both the equations 

of Reference 1 (for which B = g = 0 , etc. ) and the equations of Table 

o _ o o 

T LT G 

4 herein. The following material constants were employed. 


E f 

(Young's 

modulus of filaments) 

E b 

( Y oung's 

modulus of binder) 


(Pois son 1 

1 s ratio of filaments) 

^b 

(Pois son 1 

1 s ratio of binder) 


10 


21 


x E 


f 


0 . 2 


0. 35 


These values correspond to those for E-Glass in epoxy. Values of 3 were 

found for these material constants by derivation of the upper bound Hashin- 

Rosen (Ref. 2) elastic constants of unidirectional reinforcement at values of 

volume fraction filament of 0. 2, 0. 4, 0. 6, and 0. 8. 

Results of the calculations are given in Figure 27 for the filament orien 

tation shown. Both methods of calculation yield different compliances in the 

2-and 3-directions (values of A and A.) from that in the 1-direction. Differ 

4 o 

ences between the two methods are not substantial, as was anticipated, nor 
are the compliances found substantially at variance in the three directions. 

F or engineering purposes, the compliances in the three directions are pro- 
bably so nearly equal that differences among them can be neglected. 
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Three Phase Composite Compression Members 

In the foregoing studies of the efficiency of combination metal/composite 

structures the efficacy of 0° reinforced composite material to provide strength 
to the structure while some other material or configuration provided adequate 
transverse stiffness or continuity was evaluated. On the presumption that such 
0° reinforcement would truly be an extremely high compressive strength material, 
these evaluations were indeed favorable. Because of the relative weakness and 
lack of stiffness of available resin binder systems, however, it is not obvious 
that the 0° configuration necessarily has the maximum compressive strength. 
Rather, as this section will demonstrate, some transverse reinforcement may 
be desirable. 

Reference 4 described the mechanics of stabilization of uni -directional 
filaments in compression. As noted there, three failure modes may be encoun- 
tered: 


(1) a mode in which the shear stiffness of the binder material governs, 
described by the equation 

C - 

/ 

where 




ultimate compressive stress in filament direction 
$6 shear modulus of binder 
* volume fraction filament 

w 

(2) a mode in which the strength of the binder is critical, described by 
the equation 
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„ . /■» & A 

c Y 3(l~Vf) 


where 

<Tb„ yield stress for binder 

y 

Young's modulus for filament 

and (3) a mode in which the stiffnesses of both filaments and binder are 
operative, described by the equation 






The lowest value of found from the foregoing equations is, of course, the 

failure stress. 

In all three modes of failure the binder properties are important. Accord- 
ingly* increases in binder properties may be expected to be reflected in increases 
in compressive strength of uni -directionally reinforced composites. One method 
for increasing effective binder properties would appear to be to take some of the 
uni -directional, load -car rying filaments and distribute them uniformly through- 
out the binder. Thus, while for every failure model the value of would 

c 

tend to be reduced by the reduced volume fraction of load-carrying filament 
(v£ in above equations), the net effect should be expected to be an increase in 
due to the improvement in binder stiffnesses and strengths. 

In order to assess the potential magnitude of compressive strength 
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increase available by this approach, the assumption was made that the binder 
properties attained by the utilization of some fraction v of the filamentary 

D 

material (the original volume fraction filaments v, was assumed to be 0.5 

r o 

throughout) for binder reinforcement would be given by simple, rule -of -mixture - 
like, expressions such as 


£ b * 


where 





v 


/ 


is the Poissons ratio found by the rule of mixtures 
for filaments and binder 


and 

% * ( , -%) or 6, 

with the primes denoting the properties after binder reinforcement. 

Using the foregoing equations, three possible composite combinations were 

evaluated: (1) boron binder reinforcement for uni -dir ectionally reinforced boron/ 

epoxy; (2) glass binder reinforcement for uni -directionally reinforced glass/ 

epoxy; and (3) glass binder reinforcement for uni -dir ectionally reinforced boron/ 

epoxy. In the first two cases the total amount of reinforcing material in the 

composite was held constant, so v- + v = 0.5. In the last case v was held 

b 1 

constant and the total amount of reinforcement therefore increased as v was 

b 

inc reased. 
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The results of the calculations for these three combinations are plotted 
in Figures 28, 29, and 30. In Figure 28, curves for the boron -boron com- 


posite are presented for the following nominal constituent properties: 


v 

V 

'V 


414 GN (60,000 ksi) 



rY _ 2 
m 


3.45 

GN 

m z 

m 

(500 ksi) 

2.76 

GN 

m 2 

(400 ksi) 

0 . 1 

GN 

(15 ksi) 


0.2 
0. 35 


In Figure 29 the glass -glass combination is considered for the same binder 
properties and the following nominal glass properties 


E 


f 



72.45 GN 


2.76 GN 


0.2 


(10,500 ksi) 
(400 ksi) 


The boron-glass combination of Figure 30 uses the same nominal properties 
for each of the materials as those given above. 

For properties like those of the boron/epoxy composites of Figure 28, 
the conversion of a small percentage of the axial filaments into transverse 
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filaments which effectively stiffen the binder should be advantageous, raising 


the compressive strength from 2.75 (392 ksi) to 3.8 (550 ksi), 

m^ m^ 

approximately (see Figure 28). If the properties are more like those of the 
glass/epoxy composites of Figure 29 > however, the same gain is not to be 
expected (see Figure 29). If glass is added to boron/epoxy, on the other hand, 


as in Figure 30, substantial gains are to be anticipated. 

Because of the gains found in the calculations for Figure 28 > a series of 
uni -directional boron/epoxy compression specimens having various percentages 
of chopped -glass filaments added to the binder were fabricated and tested. The 


results are shown as the points in Figure 31, with the large and small circles 


representing nominal glass -filament lengths of 0. 16 cm. (1/16 in. ) and 0. 08 cm. 
(1/32 in. ), respectively. Despite the large scatter, strengthening appeared to 
be achieved in two cases. 


Both the scatter in the test data, and the sensitivity of strength of the 
boron/ epoxy composites of Figure 28 to the small amounts of binder reinforce- 
ment point to the critical nature of the role played by the binder in the composites 
for the development of high compressive stresses. The marked difference be- 
tween the curves of Figures 28 and 29 suggests that this criticality is substan- 
tially greater for boron than for glass reinforcement. 

To explore the sensitivity of boron composites to binder properties 


somewhat further, the calculations represented by the curves of Figure 31 were 
also made. These curves show the effects of changes in binder strengths 
and binder -filament-reinforcement strengths on the compressive 
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strength of the three -constituent composite when the "binder strength" mode is 
critical. The compressive strength is found to depend most critically on the 
binder strength at the low ratios of binder reinforcement. 
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STUDIES OF MATERIALS- PROPERTIES TESTS 


The testing of filamentary composites to determine their mechanical 
properties has proved to be more difficult than the testing of homogeneous 
materials like metals. The difficulties take several forms; they may derive 
from discontinuities encountered in specimens designed to provide a diminished 
cross-section, as is often done to avoid failures at points of load introduc- 
tion; they may be associated with the shear weakness of the binder leading 
to undesired modes of failure; or they may arise because a type of test 
known to be inadequate is used for economy or other reasons. In this section 
some of the aspects of mechanical-property testing of composites are 
examined both analytically and experimentally to help establish guidelines for 
improved techniques. 

The NOL Ring Split- Dee Tensile Test 

The N. O. L. ring "split-dee" tensile test, shown schematically in 
Figure 32, has the advantages that the specimen may be readily fabricated by 
winding, and the test may be performed in a conventional universal testing 
machine without special fixtures (other than the split dees themselves). The 
test has the disadvantages, however, that (1) no test section is available at 
which strain gages may be mounted to measure the stress- strain properties, 
and (2) more seriously, the test introduces substantial bending moments in 
the ring where the split occurs between the two dees. These bending moments 
extend above and below the split, as shown in Figure 33. Because the mem- 
brane stress in the ring causes it to increase in diameter as the load increases, 
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the ring tends to pull away from the corners of the dee, and the maximum 
bending moment is that associated with the load times the deflection to the 
center of the ring cross-section near the corners of the dees. The magnitude 
of the bending may be found from the following equations (developed in Appendix b). 



where 



may be evaluated by the trial and error solution of the 


equations 



and 




(i) 


( 2 ) 


( 3 ) 


such that 



( 4 ) 
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where 



I 

R 

o 


t 

P 

A 

1 

6 and 0 


G 


12 


bending moment at center line of ring, at split in dees 
Young's modulus in direction of filaments 
moment of inertia of ring cross-section 
initial inside radius of ring 
ring thickness 

load acting in each half of ring (one -half of applied load) 
ring cross-sectional area 
unsupported semi-span between dees 

slope at point of departure of ring from dee, measured 
relative to direction of load application, 
shear modulus along filaments 


Evaluation of equation (1) for a typical E-glass reinforced epoxy ring 

for = 7.30 cm. (2.875 in.) yields the curve of maximum bending stress 

vs. ring thickness given in Figure 34. Obviously these bending stresses are 

of sufficient magnitude to raise questions about the engineering merit of the 

p# 

tensile strength values measured as the > A stress at the maximum load on the 
ring. 

(Curiously, the magnitude of the bending moment is relatively insensi- 
tive to the composite properties. For example, the stress increment due to 
bending in an hypothetical boron/epoxy ring with E^ = 414 (60,000,000 

psi) is essentially the same as that in a glass-epoxy ring with E^ = 51.75 

m 
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(7, 500, 000 psi) at the same load. The increased bending stiffness of the boron/ 
epoxy is compensated for by the decreased diametral expansion of the ring 
under load, so that the bending stresses remain essentially constant as 
varies. ) 

Variations of the Split-Dee Test 

The high bending stresses encountered in the split-dee test are rather 
disappointing in view of the many attractive features of the test. Accordingly 
the question naturally arises as to whether minor changes may be made in 
specimen or fixture de sign which will reduce or eliminate the bending. Some 
possible variations of this nature are considered below. 

The "racetrack" specimen - The provision of a straightaway section adjacent 
to the split in the dees is a first logical step toward the improvement of the 
N. O. L. ring split-dee tensile test. Such a specimen is shown schematically 
in Figure 35. 

The analysis of the bending of the critical section of the racetrack is 
similar to that of the ring, and is described in the following equations: 
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where now 


is evaluated from 



again with 

9, = e, 

where s = length of straightaway, and the other symbols are as before. 

Results from the analysis of this racetrack specimen compared to those 
from the N. O. L. ring are compared in Figure 36. Figure 36 shows that while 
bending is not eliminated by the straightaway, it is substantially reduced, - 
even by a relatively short straight section. For example, just a 1 cm (0.4 in.) 
straightaway reduces the bending moment to less than one -half that of the 
circular ring. 
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The mechanics of the reduction of maximum bending moment are 
suggested in the (exaggerated) sketch of Figure 37. Under load the mid-point 
of the straightaway deflects inward toward the center of the track, so that 
even though the curved portions pull away from the dees the eccentricity of 
the tensile load at mid- span is small and hence the moment associated with 
it is small, - i. e. the maximum moment is still that near the corners of 
the dees, not that at mid- span. The fact that the moment in the straight 
portion is small, however, is advantageous for the use of strain gages in that 
region. 

Experimental confirmation of a qualitative nature of the foregoing 
analyses of split-dee tests were obtained by photoelastic tests (Figure 33 
and 38), and by strain measurements on an enlarged, aluminum- alloy ring 
(Figure 39). Even with the enlarged ring the stress gradients were so steep 
around the circumterence that the gages could not be located accurately enough 
to make a quantitative check of the analyses. Indeed, the strain gaging 
problems pointed up the merits of some such approach as the following for 
further improving the split- dee test. 
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Compression Tests of N.O.L. Rings 

A simple compression test for an N.O.L. ring -type specimen (analogous 
to the split-dee test for tension) would be useful for the evaluation of the com- 
pressive properties of filament -wound composites. Efforts to develop such a 
test under this contract were not successful. The specimen tried is shown 
schematically in Figure 40, consisting of two short segments cut from an N.O.L. 
ring and mounted back to back as shown. 

Application of tangential end loads to the segments of ring in Figure 40 and 
pressure -like forces from the supporting, hour -glas s -like core ideally would 


reproduce the compression induced in an entire ring by external pressure. 


That is , 


moment equilibrium at any station 



along the segment is represented 


by the equation 

PR[i-«m(<- 9,)] - f * [R*4ui{0-9,)]d8 

h J *. 

where 


P = 1/2 total load 

R = ring radius 


Of - 1/2 angle included by segment 

= intensity of normal load on ring at station 
that is 4Y X hr (&) 


load. 


As might be expected, this equation is compatible with a simple pressure 


Thus , if 


w * 
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PR I fsmO - 9 PRf ' - cos* cosO f - sinus/* Q t J 

* PR [t - ***(*”&,)] 

Tests of specimens fabricated as sketched in Figure 40, however, were not 
successful in maintaining the compatibility so readily expressible in equation 
form. Attempts were made to replace the resin filler of Figure 40 with aluminum- 
leading to a nice test of the aluminum (Fig. 41); attempts were made to test 
specimens of larger radius of curvature in this fashion, and a maximum failure 

stress of 5 9 (86 ksi) was so achieved in a glass/epoxy composite. Despite 

m^ 

the fact that failure in this case appeared to be by shearing at the quarter point 
of the specimen with no evidence of overall column instability (Fig. 42), the 
stress is less than can be achieved with well collimated, straight compression 
specimens and is probably not representative of the material strength. 

Some evidence of '’brooming" of the ends of the specimens was apparent 
even in this best test. A number of further tests of like specimens with ends 
cast in cerrobend were made to try to avoid such end failures. Stresses in no 
case were higher than the foregoing value. 

Initially straight specimens appear to be required for compression testing. 
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Low-Melting Alloy Casting Fixture of Compression Tests 


The use of a fixture in which the specimen ends can be cast in a low- 
melting point alloy has proved successful in preventing premature "brooming" 
type failure in tests of straight compression bars. Such a fixture is shown in 
Figure 42. The fixture incorporates end plates having holes, - of the same 
shape as the specimen c ros s -section , - but approximately 0.6 cm. (1/4 in.) 
oversize so that . 3 cm. of the low-melting alloy can be cast all around between 
the hole wall and the specimen. 

Conventional 3 to 1 aspect ratio, uni -dir ec tionally reinforced compres- 
sion specimens have been tested in this type of fixture with consistently high 
stresses at failure. Consistency and modes of failure both suggest that these 
high stresses are representative of the compressive properties of the material. 
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Elastic Moduli 


The desire to measure the elastic constants of a uniaxial fibrous com- 
posite by using specimens which represent the filament winding fabrication 
process has led to the use of thin-walled circumferentially wound tubes 
(e. g. Ref 29). Axial load and internal or external pressure tests of such tubes 
can be used to measure four of the five independent effective elastic con- 
stants of the composite material. (See Ref. 11 for a discussion of independent 
effective moduli of composites. ) The use of these thin cylindrical tubes to 
measure the elastic moduli of a composite laminate appears to be a most 
suitable approach. However, it is not without problems; the first of which 
is that of definition of the moduli. When a laminate is to be used as a plate 
or shell structure, the desired elastic relations may be written with respect 
to the principal geometric axes, denoted as the 1 and 2 directions, as: 



e are the three middle surface strain components 
H are the three curvature components 
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N are the three stress resultants 


M are the three moment resultants. 

In this form, for an arbitrary laminate, the C matrix is a 6 x 6 array of e- 

ij 

lastic constants which must be determined experimentally. Because the lami- 
nate is non-homogeneous there may be coupling between extension and bend- 
ing. The existance of certain non-zero terms in the C matrix complicates 
the relationships among these elastic constants and the so called engineer- 
ing constants: Young's moduli and Poisson's ratios. This is best illustrated 

by considering alternate forms of Equations (8) - (13). These Equations may 
be written as (e. g. Ref. 2): 


N. = a. . e + b H 

1 j ij j 

M. = b e + d x 
1 ij j ij J 


These Equations are represented by: 



Equation (15) can be manipulated to yield: 



and 


(14) 


(15) 


(16) 
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A 





( 17 ) 


As an example we now consider the Young's modulus in the 1 direction 
which can be defined as the ratio of average stress to strain and is therefore 
given by: 


N 


E i ■ 


t e. 


(18) 


where t is the laminate thickness. 


This modulus can be evaluated from Equation (17) for N = N = M = 0, 


2 3 


as 


E = 


1 tA 


11 


(19) 


Or it can be evaluated from Equation (16) for N = N = x = 0 as: 

2 3 i 


E 

1 



( 20 ) 


The definition of Equation (19) is consistent with the usual practice of having 
all but one of the stress components vanish. However, the restrictions assoc- 
iated with the definition of Equation (20) are representative of the loading con- 
ditions on the shell specimen. That is, X =0, because of symmetry. These 

i 

two definitions are not the same unless b . = 0; in which event, it follows that 

ij 
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f - g. . - B - C =0 and A - e 
ij iJ i j ij i j kj . 

It is suggested that the rational approach to this problem is to de- 
emphasize the calculation of Young’s moduli, etc. and to utilize a six by six 
matrix definition of the elastic constants of a laminate. When a particular 
extensional or bending stiffness is required, the definition thereof should be 
explicitly stated. 

A second aspect of the modulus problem is the influence of the form of 
the material anisotropy and the specimen configuration upon measured values. 
The filament wound cylindrical tube under consideration is a material possess 
ing cylindrical anisotropy. In the case of a circumferential or other winding 
pattern symmetric with respect to the principal geometric axes, the material 
may be considered to be cylindrically orthotropic. An analysis of the stress 
distribution in such a medium for various applied loads is available in Refer- 
ence 30, which shows that the stress distribution even for simple applied loads 

is non uniform. For example the axial stress , a , in a cylindrical tube sub- 

z 

jected to an axial load is not constant unless V = V This condition is gen 

rz Q z s 

erally not satisfied in a filament wound shell, and for a circumferentially 

wound shell, the two Poisson ratios are definitely unequal. In this case, the 

modulus defined by Equation 18 is an average value and can be related to the 

true value by using the expression for the axial stress, a . This result is 

z 

presented on page 254 of Reference 30. However, there is evidently a typo- 
graphical error in Equation (40. 4) which is not diemsnionally correct as pre- 
sented. The desired result was therefore rederived for the particular type 
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of material symmetry represented by the circumferentially wound tube. The 
extensional stresses and strains for a cylindrically orthotropic tube are 
given by: 



These relations can be further specialized for the circumferentially 
wound tube inasmuch as the material plane normal to the local fiber direction 
can be assumed to be a plane of elastic isotropy. Thus 


E 

z 


re 


= E 


r 



u 


6z 



( 23 ) 


106 



Equations (21) and (22) are also applicable to filament wound tubes which 
are symmetric with respect to the longitudinal and circumferential direc- 
tions and which have a sufficient number of laminae such that coupling between 
extension and shear is negligible. To take advantage of this wider applicabili- 
ty the analysis will not utilize the simplifications offered by Equation (23). 
Since the applied loads are symmetric, there will be no shear stresses or 
strains and the constitutive relations for this problem are those of Equation 
(21). By symmetry, all variables will be functions only of the radial coordi- 
nate. Thus the equilibrium equations simplify to: 


do a -a „ 

+ -I— §_ = o 


(24) 


which can be satisfied by selecting a stress function, cp, such that 


a 


1 dcp 
r dr 
.2 


(25) 


■V7# 

dr 

The only compatibility equation which is not satisfied automatically 

iS: d£ r d (l*H\ „ 

- r d7~ + dr r — r° 


(26) 


Substitution of Equations (21) and (25) into (26) yields a fourth order differen- 
tial equation on cp given by: 


4 IV 3 III 
r cp + 2r cp 


2 Q II 

r gcp 


2 I 

+ rp cp = 0 


(27) 
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, 2 S r 
where p = — 


<\-V V ' 
rz zr 

E 




For p / 1 the solution is : 


- - A 1-P , T, 1 + 3 , r 2^ 

- * >-r 1 t>r + Cr + D 


The stresses resulting from this are: 


(28) 


(29) 


a r = A(l-p)r' (P+1) + B(l+P)r P_1 + 2C 
a - Ap(p-i) r " (P+1) + Bp(p +l)r P _1 + 2 C 


The cylinder is considered to be subjected to an axial shortening re- 
sulting in a uniform strain, e = e, and internal and external pressures, p 
an( i P q > respectively. This train displacement relations simplify to 


e 

r 


e 

9 


_du 

dr 

u 

r 



Substitution of Equation (30) into Equation (21) defines the strains e and e 

r 

Equation (31) offer two methods of finding u. For a unique determination: 


0 * 


2C = 



P 

z 


( 32 ) 


Use of the boundary conditions: 
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( 33 ) 


a (r = r .) = - p. 
r i 1 

o (r=r ) = - p 
r o o 


defines the remaining two constants in Equation (30). This yields: 





-2p1 

N -( 3 + 1 ) 

O’ — n / \ 

3 -l 

n / r \ 

*-(t) 

1 

r P i(r.) 

_1-C 2p J 

" P o ( r ) 
o 

I 

I— 1 

1 

0 

tv 

XD 

1 


1 





1 - C 



-(3 41 ) 

r ^3-ii 

/ r \ 

1 - c 

'v 

_ 1 -C 2P _ 


r # 

where C = — — 
r 

o 


v -(3+D 

2f3 n 

1 + — 
r 

3 -l 

r x - 2 3n 

a = P.3 (^) 

0 i \ r. / 

l 

o 

- ^(f) 

O 

l 

23 

1- C J 

23 

1 - C 


P 

z 



r»c«i 

+ ./rx-« +1 » 

ri-c^i 

1 -C 2P 

+ P (— ) 
i 

i-c 2B 




l— _J 


a = E e + a + y CT 
z z 0 z 0 i*z r 


(34) 


(35) 


(36) 
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Axial Compression Test 


It is now possible to compare the results obtained by the tube tests with 
the actual elastic moduli. First, consider the tube in axial compression. The 
measured modulus is given by: 

E - ~ (37) 


Thus : 


a is the average axial stress (total load divided by 

Z ' 

tube cross-sectional area) 



a 

z 




rdr 


(38) 


Substituting Equations (36), (34) and (35) with p. = p =0, into Equation (38) 

1 o 

and the result into Equation (37) yields, after much manipulation: 


E - E 
z z 

E 


z 


(K -v ) 

6z rz 

(1 -V V ) 

rz zr 


lll± D j, 

2 2 \ 2 

(P -1) 1+P 


(hC 2g )(hC Z ) -4C 
(1-C 2B )(1-C 2 ) 



For a thin shell we find by taking the limit of Equation (39) as C-»l that 


E = E 
z z 

Thus the thin-walled circumferentially wound tube appears to be a suitable test 
technique for measuring the transverse Young's modulus. The magnitude of the 
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error associated with this technique is indicated by computing the quantity Y 
in the following equation as a function of g : 


E -E 
z z 

E 


z 


iv a -u r 

8z r z 

\-U V 
rz zr 


( 40 ) 


where the correction factor, Y, is determined by comparison of Equations (39) 
and (40). 

The factor Y was evaluated for a radius ratio of C - 0. 9 which is a rela- 
tively thick shell for this purpose. Correction factors, Y, for smaller radius 
ratios, C, will be smaller. It was found that: 


Y < 0. 01 2, < 0 £ 6 

Y < 0. 1 2 < (3 ^ 10 

It seems safe to conclude that the circumferentially wound tube in axial 
compression is a valid test for the Young's modulus transverse to the fibers of 
a uniaxial composite. 
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Internal Pressure Test 


An internal pressure test has been used (e. g. Ref. 29) to determine the 
Young's modulus in the fiber direction and the Poisson's ratio. Strains are 
measured on the inner and outer surface. For this loading we may examine 
the equations for strain by substituting Equations (34) - (36) into Equation (21). 
For zero longitudinal strain and zero external pressure this yields: 



The measured surface strains will then be given by: 



(41) 


(42) 


W = !e h 


P + 1 


1-c 


2f3 


(43) 


As in the preceding section, we choose to examine the error associated 
with the use of equations developed for isotropic materials. The alternative 
is to write the simultaneous equations relating all the measured strains to all 
the elastic moduli and solving for the latter. However, this alternative pro- 
cedure is complex and perhaps not possible as there is no apparent direct 

measurement to define the transverse shear stiffness, G . Thus, we con- 

zr 

sider 


E 


0 



(44) 
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where 


is the effective circumferential thickness for 


For equilibrium: 


E 


e 


zero axial strain. 

is average circumferential stress. 



( 45 ) 


The stiffness defined by Equation (44) can be evaluated by the use of e. g. (45) 
and either Equation (42) or (43). We denote the results Eq and E° respec- 
tively. Treating the thin shell we write: 


= 1-6 


(46) 


where 6 is the shell thickness to radius ratio. The results are: 


■ £ 4 ] 


(47) 


E n - c 


2 G rz S 0 


(48) 


Equation (47) shows that the use of the external surface strain measurement 
yields a result for the extensional stiffness which has an error measured by 
half the thickness to radius ratio, e. g. 5% error for a radius to thickness 
ratio of 10. On the other hand, the inner surface measurement involves many 
other material moduli and is therefore unsuitable for simple use in the circum- 
ferential modulus determination. 

Note also that there is a significant difference between internal and ex- 
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ternal strains. Thus: 


e s <r i ) ■ £ e (r e ) 

V r o» 



6r 


G S. 
rz 0/ 


6 

1-(P+1)6 


(49) 


Tension Test 

Reference to Section 43 of Reference 30 shows that the stress distribu- 
tion in a cylindrically anisotropic hollow shell for which there exists at each 
point a plane of elastic symmetry normal to the axis is obtained in the same 
way as in an isotropic rod. That is, a is the only non-vanishing stress 

component and 


a. 
u z 


M 

I 

o 


r 


The cir cumf erentially wound circular shell satisfies these symmetry require- 
ments and therefore provides an exact measure of the in-plane shear stiffness. 



CONCLUDING REMARKS 


The application of existing analyses of strength and stiffness of fibrous 
composites to the analysis and design of composite structures has been treat- 
ed. It has been shown that, although the present understanding of failure 
mechanisms is incomplete, it is possible to construct a rational set of fail- 
ure criteria for fibrous composite laminates. Additional studies of the fail- 
ure mechanisms of a uniaxial composite have been performed by using limit 
analysis techniques. Also the effect of stress concentrations upon tensile 
strength has been explored and initial studies were made of the complex 
moduli of fibrous composites having viscoelastic matrices. 

The design criteria were used in the evaluation of columns, plates and 
panels in which uniaxial composites for high strength were combined with 
isotropic metals or composites for high stiffness. Designs offering substantial 
improvement over those made of one material only were demonstrated. Im- 
proved composite materials were achieved through the use of a third phase as 
a local reinforcement of the matrix. Also approximately isotropic three- 
dimensional materials were designed. 

Test methods were analyzed leading to a verification of the suitability of 
thin-walled tubes for modulus measurements and the unsuitability of the NOL 
ring for tensile strength measurements. Improvements in methods for measure- 
ment of tensile and compressive strength were demonstrated. 

The studies emphasize the feasibility of using present methods for the 
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mechanical analysis of composites in preliminary design studies. They also 
indicate that our understanding of composite failure mechanisms is a continual- 
ly evolving one requiring additional theoretical and experimental study. 
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APPENDIX A - Limit Analysis For Composites With Elastic -Brittle Fibers 


Five types of surface loadings are considered: 

1. Tractions equivalent to a uniform shear stress T are applied on 
the entire boundary surface. 

Use will be made of the theorems of limit analysis which will not be 
stated here. For lower bound construction, a uniform stress field 

0 TO 

o 

T 0 0 

o 

0 0 0 

is chosen as a statically admissible stress field where T is such that Von 

o 

Mises 'yield criterion of the matrix is not violated. Then it can easily be 
shown that 



T.. = 
U 


where is the yield stress of the matrix in simple shear. 

For upper bound construction, a kinematically admissible velocity field 
is chosen as follows: 

(a) In the region of the composite specimen not occupied by composite 
cylinders and on the boundary of the composite cylinders. 


“i =0 ’ 


U 2 = Y 0=Y 



(2) 
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where v is 
0 


any real number. 


(b) In any composite cylinder, the velocity field ^ is the elastic dis- 
placement solution to the displacement boundary value problem with displace- 
ment boundary conditions (Eq. 2) prescribed as formulated in Appendix 2 of 


(Ref. 11) with the modification that the fibrous core is rigid. 

An application of the upper bound theorem gives an upper bound (t L ) 

12 u 

as a function of p and the latter being the volume fraction of all the composite 
cylinders embedded in the composite specimen: 


(^1 = V Vi “f *’ 

where 


(3) 


I = 
1 


-^jT 

nd-P ) J p J o 


+ 2 


cos 0 d0dR 


R 


R 


In the case of ''random array" (Ref. 11) for which v = 0 and the fiber 

Z 

volume fraction of the composite specimen = p , Equation (3) becomes 



2. Tractions equivalent to a uniform shear stress are applied on 
the surface of the specimen. 

The condition of plane strain is assumed. Then, the Von Mises 1 yield 
criterion for the matrix reduces to a simple form 
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( T Z2 ' T 33f + 4T 23 S 4k b 
For lower bound construction, a uniform stress field 
0 0 0 “ 

0 0 T 

o 

0 T 0 

o 

is chosen as a statically admissible stress field where T is such that Eq. (5) 
is nowhere violated. Then it follows that the lower bound for the limit load 



which is independent of fiber volume fraction. 

For upper bound construction, the same principle is used here. A 
kinematically admissible velocity field is chosen as follows; 

(a) In the region of the composite specimen not occupied by composite 
cylinders and on the boundary of the composite cylinders, 



u 


1 


= 0 , 



u 

3 



(6) 


(b) In any composite cylinder, the velocity field u is the elastic dis- 
placement solution to the displacement boundary value problem with displace- 
ment boundary conditions (Eq. (6) ) prescribed as formulated in Appendix 1 
of Ref. 11 with an additional condition that the fibrous core is rigid and the 
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binder shell is incompressible. 


An application of the upper bound theorem gives an upper bound 


as a function of 3 and v : 

1 





k bYV» 


(7) 


where 


fa C 


J 2 2. 3 


"( 1-0 ) •'3 


R) + ♦ (p, R) cos 0 ded R 


and 


ll / (B. R) = — 
1 R 3 4 


z 

(e\e 2 + ,)-l(e 2 +i)-Cl + [.ilMllL + 3f ? R z 

“ “ ' r‘ j L 2 




R L 


(0 4 +0 2 +i) _!_( 3 2 +1 )!^ 


R 


(4B 4 +3 2 +i) 


2 2 
+ 33 R 


3. Tractions equivalent to uniform axial tension x on the composite 
specimen under the conditions of plane strain. 

Using the same principle, the lower bound for the limit load is 



For upper bound construction, a kinematically admissible velocity field 
is chosen which is obtained from the one constructed for the preceeding case 
through an orthogonal transformation such that 
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ooo 


u = d u 1 with d = 0 — — 

1 ij J ij /Z fz 


0 




where u! are the velocity components used in the preceeding case referred to 
an x' -system. Then, after some manipulation, it turns out that 



Another kinematically admissible velocity field was constructed: 

I ft the region of the composite specimen not occupied by composite 
cylinders and on the boundary of the composite cylinders, 



or in polar coordinates, 

Y o Y 

U = o, u = r cos 29, u = - r sin20 
1 r 2 b 6 2b 


b. In any composite cylinder, the velocity field in nolar coordinates 


is assumed to have the following form: 


V 

u = — ^ - p(r) cos 2 9 
r 2 

for r - r £ r 

Y f b 

Ug = - Q (r) sin 2 0 
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where p (r) and Q(r) are arbitrary functions and 


p(r b ) = Q ( r b ) = r b 


( 8 ) 


In addition, we assume that p(r^) - fi(r^) = 0. Becuase of the assumed 
incompressibility of velocity field, p(r) and fl(r) satisfy the following equation: 


dp(r) 

dr 


~> n I 

+ kv g / _ q 


(9) 


2 3 

Furthermore, let p (r) = A^(r-r^) + A^(r-r^) 


A = 
1 


r + 2 r, 
__f b 

(r - r 
b f 


where the constants 


A = - 
2 


V r b 

<r b-V : 


have been determined from (8) and (9). 

Now the velocity field of the entire specimen is constructed and the 
upper bound theorem is applied to obtain a new (j | i n terms of fiber- 

\ zz Ju 

volume fraction for the case of "random array". Incidentally, this same vel- 
ocity field, after an orthogonal transformation, can also be used to obtain a 

'U 

volume fractions between 0 and 1, both 


new^T 1 . Unfortunately, numerical calculations show that for all fiber- 


( t 4u 


and (t ] obtained here are 
23 L 


slightly higher than their corresponding values obtained previously. 
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4. Tractions equivalent to biaxial uniform tension t and t 

22 33 

(t ^ T 33 ^ on t ^ ie composite specimen. 

For definiteness assume > Then following the ideas used in 

Case 3, it can be shown that 




5. The type of surface loading considered here is for T in the follow- 


ing form: 


T = 

ij 


r o 

12 


T T 

12 22 


( 10 ) 


where the constant stress components r and r are related in the following 

12 22 

way: 


T 


12 


a r 


22 


(ID 


with a ^ 0. This amounts to a proportional loading of combined uniform shear 
stresses r j2 and uniform tensile stresses t on the boundary surface S of 
the specimen. 

Since both T ^ and T are assumed finite, it is obvious that a= 0 
corresponds to the case where only uniform tensile stresses are present. 
On the other hand, a-» ® corresponds to the case where the specimen is sub- 
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jected to uniform shear stresses r 

12 

In lower bound construction, a uniform stress field of the form (2) 
supplemented by (3) is chosen as a statically admissible stress field to obtain 
lower bounds for the limit loads and > respectively. It turns out that, 


and 



for 0. 


( 12 ) 


IT rom Equation (12) it is obvious that for the special case where a - 0 


(i. e. for uniform transverse stress r acting alone), 



(13) 


In the case where a-+ 00 , corresponding to the case where the specimen 


is subjected to the uniform in-plane shear stresses t alone, 

12 


lim ( T 12 ); 
a -» “ k^ 


= 1 
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For upper bound construction, a kinematically admissible velocity 

field u is chosen to be a linear combination of the two kinematically ad- 

-4 

mis sible velocity fields used to obtai n(r^\ and(VV) > respectively: 

v 12 'u 22 u 

(a) In the region of the composite specimen not occupied by composite 
cylinders and on the boundary of the composite cylinders: 


( 1 ) ^ ( 2 ) 
u = u + u 


(14) 


where 


and 


= (0, v x , 0) 
-► 11 


(2) (n 

u = ^0. 


2 

T 


2 ’ 


2 X 3; 


y = uu Y 
1 2 


(15) 


where (JU is any real number. 

(b) In any composite cylinder, the velocity field is the elastic displace- 
ment solution to the displacement boundary conditions ( 14 together with (15) 
prescribed with an additional condition that the fibrous core is rigid and the 
binder material is incompressible. In fact, the solution to this displacement 
boundary value problem can be obtained from the associated solutions to the 
individual problems connected with obtaining and by the prin- 

ciple of superposition. 

Then an application of the upper bound theorem gives, for the case of 
constant j3 throughout the specimen 
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( 16 ) 



where 


I 

3 


.. -2. 3 
TTU-P ) 



n (UJ, 3, R) + 0 (ua, 3, R) cos 0 + Q (3, R) cos 20d0dR 

VC, J 


G (io, 3. R) = 4 



+ 


2 2 
3 3 R 


r 4 2 

V43 40 + 1 



+ 




4 “ 

n_ ( w, 3. R) = 2 UU 2 ( 1-3 2 ) -£-=■ 

2 R 


and 


fi (3, R) = 8 £- 
J R 


r 2 -i 

LA 2 \ 3 . d 2 3 

(p 43 4lj - — P +1 ^2 


3 P 2 r 2 - lie ie ±» 


4 - 2 -] 


In the case of "random array", (16) reduces to 


/ L\ 

^ 22 Aj = 

b -| + auu 


(17) 
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Since a J in (15) is arbitrary, the lowest upper bound among the class 
of upper bounds in (17) will be obtained by minimizing the right-hand side of 
(17) withrespect to uu . 

Thus, 



min 

M 


k I 
b 3 


1 

- +auu 


( 18 ) 


and 



k a I 
b 3 

i +a ” 


will be chosen as the upper bounds for 7 ^ and r respectively. 

22 12 


Numerical calculation is performed to obtain 



and 



from 


(18) for different values of 3 and a. It is interesting to note that in the num- 
erical calculation, for any giver UL) which minimizes the right hand side of 
Equation (18) is a monotonic increasing function of a but (ju ^ a (except when 
a - 0, then a) = a = CM . 
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APPENDIX B - Derivation of Equations for N.O.L. Ring Split- Dee Test Analysis 

The assumptions are made that the N. O. L. ring exactly fits around 
the two dees at zero load and that there is no friction between dees and ring, so 
that the portion of the ring that remains in contact with the dees is subjected 
to essentially a radial pressure which induces a hoop tensile force P in the 
ring. This hoop tension produces a circumferential expansion ^ of the ring, 
equivalent to: 

A c = — fer)(R e + Y) ,i> 

where 

E^ Young's modulus in hoop direction of ring 

A cross-sectional area of ring 

initial inside radius of ring 
t r ing thickness 

This circumferential expansion is converted to a straight-line separation 
between the dees, so that 

very nearly, and the problem is that of the analysis of a beam at the split in the 
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dees as sketched below: 



X 


That is, the beam is tangent to the dee at some point at a distance 

from the center line of the split between the Dees. At the centerline the beam 

is subjected to the tensile force P and a moment . At the point of tangency 

the force P is also acting and a fictitious moment M is hypothesized of the 

magnitude required to form the ring to its initial radius of curvature R + — . 

o 2 

So 





R *•£ 


In other words the analysis considers the beam as initially straight when un- 
stressed. The desired bending moment at the center line will thus be found 

by subtraction of the fictitious moment M from . 

o L, 

Selecting x-y coordinates as shown, we note that the moment M at any 
section of the beam is 


M = -Pfi-y) + M 0 


1 32 


or 


Let 



By substitution (4) becomes 


c/*y 



-k l S + 



The solution of (6) is 



( 4 ) 


(5) 


( 6 ) 


(7> 


Accordingly the slope ©. at the point of tangency of the beam to the de 


es is 


The deflection & is similarly found as 


S'* f (*-&)$ 


( 8 ) 


<r 




(9) 


133 



Combining (8) and (9) 




( 10 ) 


Because the beam is tangent to the dee at x = Jl , the angle e t subtended by 
the arc of the dee between the split and the point of contact may be found by 


trigonometry as 









( 11 ) 


And 



( 12 ) 


Finally, from (3) for y = 0 

M = M. = M 0 -PS 

Or by substitution of (9) 


(13) 



( 14 ) 
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Shear deflections may be included in the analysis by re-writing equation (7) as 



For a rectangular section n - 1.2, and (15) becomes 



Analysis to include the effect of the straightaway on the racetrack specimen 

follows the same form as above, leading directly to the equations given in the 
text. 
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/- Direction Filaments 



A7 - GbP'A* V f]+ Yf] 

A 4 1 G b[>-P.y f ]* 


Table 1. - Generalized equations for compliances of composites having uni- 
directional reinforcing filaments in the 1 -direction (from Ref. 1). 
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Table 2. 


h Direction Filaments 

Interrelationships among constants /hr 
transverse isotropy 



A g = A 4 

A * A = 

4 6 

Af “ Ag * G/g 




- Relationships among elastic constants and compliances for uni-direc- 
tional reinforcement for use in evaluations of the various j$ T s in 
Table 1 (from Ref. 1). 
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Table 3. 



- Equations for the transverse effectiveness factors 0) in terms of 
the elastic constants for uni-directional reinforcement (equations of 
Table Z solved in terms of g for Q m * values equal to unity). 
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T a ble 4. - Generalized equations for the compliances of three -dimensionally 
reinforced composites having three orthogonal planes of symmetry. 
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, [ * fEf t vEi I . 

[ %E fi _ ^Ej, L v 

I £/ / ) 

j^fej . 

^ j#//*y cm*ji ♦ cottyjcosy co/a * ^^/\ eft, v| •u 




(/+*iX*-2n) 






>y w*fl*cot i ^jco/^ co/IL t cot 


Table 4. - (Cont. ) 
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T able 4. - (Cont. ) 
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Table 4. - (Cont. ) 
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STRESS, KSI 



Figure 1. Calculated Stress -strain Curves for E -Glass and Epoxy Composite 
Laminates. 







Figure 3. Composite Specimen for Limit Analysis. (Fibers in 


X j -Direction). 
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47 


igure 4. Bounds on the Limit Load for In-plane Shear. 



. FOR ALL )9 2 (0<£ 2 <l ) 
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Figure 6. Bounds on the Interaction Curves for In-plane Shear and Transverse Tension for Various 
Fiber Volume Fractions. 
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igure 11. Expected Number of Groups of Fractures as a Function of Applied Load (Ref. 
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xpected Number of Groups of Fractures as a Function of Applied Load (Ref. 
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Observed and Predicted Multiple -Fracture Groups (Ref. 
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Observed and Predicted Multiple -Fracture Groups (Ref. 
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igure 15. Static and Dynamic Values of P , and P , as a Function of Applied Load (Ref. 4, Series B). 
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ected Number of Fracture Groups as a Function of Applied Load (Ref. 17 Tests). 
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igure 17. Variation of E with Applied Stress for Various Composite Lengths (Ref. 4, Series B). 









Scalloped Triangle Tube 





Figure 22. Changes in Column -Binding Characteristics of Scalloped Thin 
Walled Tubes at Constant R / 1 with Angle Included by Scallop, 
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Figure 23. Efficiencies of 7075-T6 Aluminum -Alloy Scalloped Triangle Thin-Walled Tube Columns 
Reinforced at Apexes with 0° Boron/Epoxy Stiffeners, and Comparison with Minimum- 
Weight Reinforced 7075-T6 Round Tubes. 
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Figure 24, Efficiencies of Beryllium Scalloped Triangle Thin-Walled Tube Columns Reinforced at 

Apexes with 0 Boron/Epoxy Stiffeners, and Comparison with Minimum- Weight Reinforced 
Beryllium Round Tubes, 




167 


Figure 25. Efficiencies of Solid-Sandwich Plates having 0 Reinforced 
of Various Materials. 
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Figure 26. Comparative Efficiencies of Wing Box Beams having ±45 Boron Reinforced Epoxy Compres- 
sion Skins and Beryllium Skins on 0° Reinforced Boron/Epoxy Z -Section Stiffener s . 


Ca/tropic Reinforcement 
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Calculated Compliances A ^ , and A^, Caltropic Reinforcement 

in the Directions of the Orthotropic axis 1-, 2-, 3- Respectively, for 
Equal and Unequal Assumed Transverse Effectivenesses (}. (E^. = 10, 

E b = n- v t= °- 2 ' ■' b = 0.35). 
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Figure 27. 




Figure 28. Calculated Compressive Strengths for Three Indicated Failure Modes 
tor a Boron/Epoxy Composite of Nominal Constituent Properties with 
Constant Total Reinforcement Volume Fraction of 50% but Varying 
Proportions of the Reinforcement Uni -dir ectional and Randomly Dis- 
persed in the Binder. 
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Calculated Compressive Strengths for Three Indicated Failure Modes 
for Glass /Epoxy Composites of Nominal Constituent Properties with 
Constant Total Reinforcement Volume Fraction of 50%, but Varying 
Proportions of the Reinforcement Uni-directional and Randomly Dis- 
persed in the Binder. 


Figure 29. 
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Calculated Compressive Strengths for Glas s -Boron/ Epoxy Com- 
posites of Nominal Constituent Properties with Constant Uni-Direc 
tional Boron Reinforcement Volume Fraction Vj = 0. 5 and Varying 
Quantities of Randomly Dispersed Glass Filaments in the Binder. 



Figure 31. Experimental Results for Tests of Glass -Boron/ Epoxy Composites 

having Uni -dir ectional Boron Filamentary Reinforcement of Nominally 
50% by Volume, and Comparison with Calculations for Various 
Glas s -Filament and Binder Strengths. 
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Figure 32. 


Schematic of N. O. L. Ring, "Split-Dee' 


Tensile Test. 
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Figure 33. Photoelastic Study of Stresses in the Vicinity of the Split between the Dees in an N. O. 
Ring "Split-Dee" Tensile Test. 






Figure 35. "Race-Track", Filament- Wound Tensile Specimen. 


177 





Ratios of Maximum Bending Moments Induced in Race-Track 
and Circular 0. 15 cm^jO, 06 in. ) Thick Split-Dee Specimens of 

(100 ksi) Axial Stress. 


Glass/Epoxy at 0. 69 


GN 
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Figure 36. 
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Figure 37. Schematic Representative of Mechanics of Deflection of "Race- 
Track" Specimen. 
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Figure 38. Photoelastic Study of Stresses in the Straightaway of a "Race- 
Track" Specimen. 
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Figure 39. Strain Measurements Near the Split in the Dee in N.O. L. Ring 
Type Split-Dee Tension Test. 
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Figure 40. Proposed Compression Specimen Made from Segments of NOL 
Ring. 
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Aluminum Insert from Specimen Simila 
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figure 42. N. O. L. Ring Segment Compression Test Specimen Showing Shear Failure. 



Figure 43. Top and End Views of Compression Specimen Cast in Cerrobend 
End Fixtures. 
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